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INTRODUCTION 


This  report  covers  millimeter  and  microwave  dielectric  measurement 
work  performed  at  the  Naval  Weapons  Center  (NWC).  The  apparatus  used 
for  obtaining  elevated  temperatures  and  making  the  measurements  is 
described,  and  the  principles  of  operation  are  set  forth  along  with  the 
results  for  various  materials. 

Missiles  are  being  required  to  fly  at  ever  higher  mach  numbers  in 
more  severe  flight  profiles  with  consequent  aerodynamic  heating  and 
potential  flight  erosion.  To  design  and  build  a  usable  radome ,  it  is 
necessary  to  know  the  propert ies--e 1 ect r i ca  1 ,  mechanical,  and  ther¬ 
mal — of  materials  used  (or  potentially  useful)  in  radome  construction 
at  the  temperatures  encountered. 

It  is  also  expected  that  missile  systems  will  operate  at  milli¬ 
meter  as  well  as  at  microwave  wavelengths.  Because  of  the  wide  variety 
of  materials  used  (or  potentially  useful)  for  radome  construction,  at 
least  some  can  be  expected  to  show  significant  change  in  their  dielec¬ 
tric  properties  over  the  millimeter-microwave  frequency  ranges. 

The  work  reported  here  represents  an  effort  to  supply  laboratory 
measurements  of  the  complex  dieletric  constants  of  materials  over  both 
millimeter  and  microwave  frequencies  at  temperatures  that  might  be 
encountered  in  a  missile  flight. 

No  single  measurement  method  appears  suitable  for  the  complete 
range — ceramics,  plastics,  laminates,  fiber-reinforced  teflons  and 
epoxies,  etc. — of  proposed  materials  over  such  a  wide  frequency  range. 
Differences  in  dielectric  constant  magnitude,  loss  tangents,  coeffi¬ 
cients  of  thermal  expansion,  machi nabi 1 i ty ,  availability,  chemical 
decomposition  rates,  and  end  products  at  high  temperatures  seem  to 
preclude  this.  In  any  case,  microwave  dielectric  constant  and  loss 
tangent  measurements  are  difficult  enough,  especially  at  high  tempera¬ 
tures,  with  enough  sources  of  error  and  uncertainties  that  it  is  impor¬ 
tant  to  have  overlapping  measurement  methods  to  supply  corroboration  or 
to  pinpoint  discrepancies. 


NWC  TP  6643 


Millimeter  and  microwave  dielectric  properties  measurement  tech¬ 
niques  can  be  divided  into  groups  (Reference  1-1).  Transmission  line 
methods  look  at  a  wave  incident  on  a  sample,  and  the  interaction  of  the 
wave  with  the  sample  is  determined.  For  example,  the  complex  reflec¬ 
tion  or  transmission  coefficient  can  be  measured,  or,  as  in  the  Massa¬ 
chusetts  Institute  of  Technology  (MIT)  Radiation  Laboratory  measure¬ 
ments,  the  reflection  from  a  sample  backed  by  a  short  can  be  determined 
(References  1-2  through  1-4).  If  swept  frequency  measurements  are 
made,  the  frequency  separation  between  transmission  maxima  can  be  used 
to  determine  the  real  part  of  the  dielectric  constant  since  these 
maxima  occur  when  the  sample  thickness  is  an  integral  number  of  half 
wavelengths  in  the  material. 

Geometries  which  employ  a  washer-shaped  sample  in  a  coaxial  line, 
a  slab  across  a  waveguide  cross  section,  or  a  flat  panel  located 
between  waveguide  horns  (see  Section  3)  are  all  examples  of  transmis¬ 
sion  line  methods  (References  1-6,  1-6).  They  have  the  advantage  of 
conceptual  simplicity  and  seeming  ease  of  implementation.  They  gener¬ 
ally,  however,  have  some  problems  in  determining  loss  tangent  in  low- 
loss  materials,  since  the  wave,  neglecting  multiple  reflections  within 
the  sample  (which  rapidly  die  off  for  the  relatively  low  dielectric 
constants  of  radome  materials),  basically  passes  through  the  sample 
only  once  or  twice.  Therefore,  the  sample  absorbs  only  a  small  frac¬ 
tion  of  the  available  power,  and  precise  measurement  of  this  quantity 
can  be  difficult.  When  making  transmission  measurements,  as  in  the 
case  of  classical  measurements  at  Brewsters  angle,  the  wave  passes 
through  the  sample  only  once  since  the  boundary  reflection  coefficient 
is  zero;  unless  very  thick  samples  are  available,  this  represents  the 
least  sensitivity  to  loss  tangent. 

These  measurement  problems  involving  sample  losses  can  be  allevi¬ 
ated  where  practicable  by  using  sources  which  are  stable  with  respect 
to  frequency  and  power  output  or  very  thick  samples,  together  with  a 
careful  accounting  for  extraneous  effects  such  as  residual  reflections 
within  the  components  of  the  measuring  system  and  (for  free-space 
bridges)  wave-front  curvature.  It  probably  remains  true  that 

transmission-line-type  measurements  are  best  suited  for  medium-  rather 
than  low-loss  material  loss  tangent  determination.  In  the  case  of 
radome  materials,  it  is  also  true  that  loss  tangents  up  to  0.01  are 
generally  considered  acceptable  and  the  precise  value  of  much  lower 
loss  tangents  is  relatively  unimportant  since  the  effect  on  electro¬ 
magnetic  performance  will  be  small. 

Resonant  cavity  methods  enclose  the  dielectric  sample  in  a  micro- 
wave  resonator,  filling  some  fraction  (up  to  100%)  of  the  resonator 
volume  (References  1-7,  1-8).  The  quantities  measured  are  usually  the 

resonant  frequency  or  the  change  in  resonant  frequency  and  the  cavity  Q 
with  the  cavity  empty  and  with  the  sample  in  place.  Another  possibil¬ 
ity  lies  in  slightly  changing  a  resonator  dimension  with  insertion  of 
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the  sample  to  maintain  the  resonant  frequency  constant.  The  dimension 
change  can  be  used  to  determine  dielectric  constant. 

There  can  also  be  problems  with  resonant  cavity  measurements.  The 
sample  filling  most  or  all  of  the  cavity  volume  is  roughly  equivalent 
to  multiplying  the  cavity  dimensions  by  the  index  of  refraction.  This 
will  greatly  increase  the  density  of  excitable  modes  for  the  higher 
dielectric  constant  materials.  Of  course,  it  is  hoped  that  symmetry 
will  rule  out  many  of  the  unwanted  modes;  however,  any  inhomogeneity  in 
sample  dielectric  properties  can  destroy  this  symmetry  and  bring  about 
the  unwanted  coupling.  Many  radome  materials,  such  as  ceramics,  com¬ 
posites,  laminates,  etc.,  turn  out  to  be  somewhat  inhomogeneous  because 
of  variations  in  the  manufacturing  processes. 

Another  problem  lies  in  the  relation  between  the  loss  tangent  and 
the  cavity  Q.  For  other  than  very  low-loss  materials,  the  cavity  Q 
will  be  dominated  by  the  dielectric  Q  (inverse  loss  tangent)  and  many 
radome  materials  have  sufficient  loss  to  reduce  the  measured  Q  below 
acceptable  limits  for  accurate  measurements.  One  can,  as  we  did,  only 
partially  fill  the  resonator  with  dielectric.  This  raises  the  Q  and 
the  measurement  accuracy  but  can  make  it  difficult  to  relate  the  mea¬ 
sured  Qs  to  the  actual  dielectric  loss  tangent.  This  is  basically 
because  the  field  distribution  in  a  partially  dielectric-filled  cavity 
can  be  very  different  from  that  in  the  empty  resonator.  So  far,  this 
is  a  problem  we  have  not  solved  satisfactorily  in  our  X-band  partially 
filled  rectangular  cavity  measurements. 

Lumped  element  measurements  form  another  measurement  category 
(Reference  1-9).  Here  a  sample  small  in  dimension  with  respect  to  a 
wavelength  is  considered  (perhaps  with  corrections)  as  a  static  capaci¬ 
tor.  For  example,  a  small  dielectric  disk  terminating  a  coaxial  line 
(Reference  1-10)  or  positioned  in  the  gap  in  a  re-entrant  coaxial  cavi¬ 
ty  (Reference  1-11)  can  be  treated  in  this  manner.  Because  of  the 
restrictions  on  sample  size  and  wavelength,  such  measurements  are  dif¬ 
ficult  at  the  higher  frequencies  and  have  been  confined  mostly  to 
S-band  and  below. 

Perturbation  measurements  are  also  employed  (Reference  1-12). 
Here  a  small  sample  is  inserted  into  a  resonant  cavity,  and  the  shift 
in  resonant  frequency  and  change  in  Q  are  measured  and  related  to  the 
complex  dielectric  constant.  The  theory  is  developed  on  the  premise 
that  the  electromagnetic  (KM)  fields  are  unchanged  in  form  and  only 
slightly  perturbed.  Such  measurements  are  attractive  when  only  small 
samples  are  available.  They  can  also  be  attractive  from  the  standpoint 
that  in  elevated  temperature  measurements  only  a  small  sample  must  be 
heated.  As  usual,  there  are  some  basic  pitfalls  and  difficulties.  If 
the  sample  is  small  enough  to  leave  the  form  of  the  KM  fields  unal¬ 
tered,  then  the  shifts  in  resonant  frequency  and  Q  are  small,  which 
leads  to  measurement  inaccuracies.  On  the  other  hand,  if  the 
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perturbation  is  large  enough  to  give  useful  values  of  these  shifts,  the 
validity  of  the  theory  can  come  into  doubt.  It  is  also  true  that  as 
the  sample  is  made  smaller,  the  tolerances  on  the  accuracy  of  the 
sample  dimensions  become  increasingly  important  and  small  dimensional 
errors  can  cause  large  errors  in  the  measured  complex  permittivity. 

Another  problem  can  lie  in  the  depolarization  factor  (Reference 
1-13).  Dielectric  bodies  in  an  external  electric  field  set  up  internal 
dipole  fields  in  opposition  to  the  external  field.  In  the  case  of  a 
small  sample  in  a  microwave  cavity,  this  will  reduce  the  perturbation 

effects  and  hence  the  resonant  frequency  shift  and  change  in  Q.  The 

depolarization  effect  is  sample-shape  and  elect ric-f ie ld-di rect ion 
dependent.  Very  long,  thin  needles  with  the  field  parallel  to  the  long 
axis  have  depolarization  factors  approaching  zero.  Flat  disks  with  the 
field  perpendicular  to  the  plane  of  the  disk  have  factors  approaching 
1.  Spheres  have  a  depolarization  factor  of  1/3  [all  in  meter-kilogram- 
second  (MKS)  units]. 

The  higher  the  dielectric  constant  of  the  sample  being  measured, 
the  smaller  the  relative  dielectric  constant  (apparent  dielectric  con¬ 
stant  divided  by  the  actual  dielectric  constant)  and  the  less  accurate 
the  measurement.  This  problem  can  be  at  least  partially  overcome  by  a 

careful  choice  of  sample  geometry,  cavity  configuration  and  mode,  and 

sample  placement.  The  depolarization  can  also  be  accounted  for  in  the 
perturbation  mathematics,  although  the  frequency  measurement  accuracy 
may  not  be  sufficient  to  take  advantage  of  this  approach. 

Cavity  perturbation  measurements  are  perhaps  at  their  best  in 
relative  measurements  (e.g.,  determining  the  complex  dielectric  con¬ 
stant  at  an  elevated  temperature  relative  to  its  room-temperature 
value)  rather  than  the  determination  of  absolute  value. 

This  report  covers  our  efforts  to  date  using  resonant  cavities 
(the  partially  dielectric-filled,  half-height  rectangular  X-band  cavity 
and  the  35-gigahertz  Fabry-Perot  open  resonator)  and  the  transmission 
line  (flat  panel  between  35-gigahertz  focusing  horns)  methods.  This 
work  is  far  from  complete,  and  it  is  our  intention  and  hope  to  pursue 
this  effort  until  we  can  perform  accurate  millimeter  and  microwave 
material  properties  measurements  on  a  wide  range  of  materials  over  a 
sizeable  temperature  range. 


NWC  TP  6643 


REFERENCES 

1-1.  "High  Frequency  Dielectric  Measurement,"  in  Proceedings  of  a 
Tutorial  Conference  on  Measurement  of  High  Frequency  Dielectric 
Properties  of  Materials,  National  Physical  Laboratory,  Tedding- 
ton,  U.K.,  27-29  March  1972.  J.  Chamberlain  and  G.  W.  Chantry, 
eds.  Place,  IPC  Science  and  Technology  Press,  Ltd.,  1973. 

1-2.  B.  A.  Burdick,  T.  J.  Lyon,  and  J.  E.  Pippin.  "Measurements  of 
Large  Dielectric  Constants  and  Loss  Tangents  at  35  Gc/s,”  IEEE 
Trans.  Instrum.  Control,  Vol.  IM-13  (1964),  pp.  318-23. 

1-3.  E.  S.  Sabisky  and  H.  J.  Gerritsen.  "Measurement  of  the  Dielec¬ 
tric  Constant  of  Rutile  (Ti02)  at  Microwave  Frequencies  Between 
4.2°  and  300°K,"  J.  Appl.  Phys.,  Vol.  33  (1962),  pp.  1450-53. 

1-4.  University  of  Michigan  Research  Institute.  A  Graphical  Method 
for  Measuring  Dielectric  Constants  at  Microwave  Frequencies,  by 
C.  B.  Sharpe.  Ann  Arbor,  Mich.,  Electronic  Defense  Group, 
Department  of  Electrical  Engineering,  1959.  (Technical  Report 
No.  94.) 

1-5.  K.  H.  Breeden,  J.  B.  Langley,  and  A.  P.  Sheppard.  "Complex  Per¬ 
mittivity  Measurements  at  Millimeter  Wavelengths,"  in  Dielectric 
Materials,  Measurements,  and  Applications  Conference  Publication 
No.  67,  IEE  London,  July  1970. 

1-6.  I.  Ye  Arsayev.  "Effect  of  the  Proximity  of  Receiving  and  Trans¬ 
mitting  Antennas  on  the  Measurement  of  Dielectric  Constant  of  a 
Substance,"  Radio  Eng.  Electron.  Phys.,  Vol.  15  (1970),  pp. 

1806-14. 

1-7.  Dielectric  Materials  and  Applications,  A.  R.  Von  Hippel,  ed.  New 
York,  John  Wiley  and  Sons,  Inc.,  1954.  Pp.  70-71. 

1-8.  R.  N.  Clark  and  C.  B.  Rosenberg.  "Fabry-Perot  and  Open  Resona¬ 
tors  at  Microwave  and  Millimetre  Wave  Frequencies,  2-300  GHz," 
J.  Phys.  E:  Scl.  Instrum.,  Vol.  15  (1982),  pp.  9-24. 

1-9.  Dielectric  Materials  and  Applications,  A.  R.  Von  Hippel,  ed.  New 
York,  John  Wiley  and  Sons,  Inc..,  1954.  Pp.  47-62. 


1-10.  Dielectric  Materials  and  Applications,  A.  R.  Von  Hippel,  ed.  New 
York,  John  Wiley  and  Sons,  Inc.,  1954.  P.  69. 

1-11.  A.  Kaczkowski  and  A.  Milewski.  "High-Accuracy  Wide-Range  Mea¬ 
surement  Method  for  Determination  of  Complex  Permittivity  in  Re¬ 
entrant  Cavity:  Part  A  -  Theoretical  Analysis  of  the  Method," 

Microwave  Theory  Tech.,  Vol.  MTT-28  (1980),  pp.  225-28;  Part  B  - 
Experimental  Analysis  of  Measurement  Errors,"  ibid . ,  pp.  228-31. 

1-12.  G.  Birnbaum  and  J.  Franeau.  "Measurement  of  the  Dielectric  Con¬ 
stant  and  Loss  of  Solids  and  Liquids  by  a  Cavity  Perturbation 
Method,"  J.  Appl.  Phys.,  Vol.  20  (1949),  p.  817. 

1-13.  C.  Kittel.  Introduction  to  Solid  State  Physics,  3rd  ed.  New 
York,  John  Wiley  and  Sons,  Inc.,  1966.  P.  378. 


NWC  TP  6643 


2 

THE  PARTIALLY  DIELECTRIC-FILLED  X-BAND  RECTANGULAR  CAVITY 


This  measurement  approach  employs  an  X-band  rectangular  transmis¬ 
sion  cavity,  as  shown  in  Figure  2-1.  Cavity  coupling  was  through 

circular  irises  at  each  end,  as  shown  in  Figure  2-2.  A  thermocouple 
embedded  in  the  cavity  wall  provides  an  input  to  a  temperature 
controller  that  adjusts  the  oven  current  to  maintain  the  desired 

oven-sample  temperature. 

The  cavity  is  slid  into  a  square  (X-band  flange  size)  channel  from 
the  ends  of  the  cylindrical  oven,  and  long  stainless-steel  waveguide 
sections  are  used  for  thermal  isolation  on  either  side  of  the  cavity. 
This  oven  is  capable  of  increasing  the  cavity  and  sample  temperatures 
to  800°F.  In  order  to  avoid  long  line  effects  from  the  inevitable 
minor  reflections  from  the  waveguide-coaxial  transitions,  etc.,  a 

rotary  vane  attenuator  was  placed  at  the  stainless-steel  waveguide 
inputs  to  provide  isolation  or  padding. 

A  receiver  (either  a  Pacific  Measurements  scalar  analyzer  or  a 

Scientific  Atlanta  receiver  with  a  chart  recorder)  was  used  to  deter¬ 
mine  the  cavity  output  magnitude  versus  frequency.  This  shows  the 
resonances  with  their  frequencies  and  Qs. 

Our  initial  efforts  employed  a  simple  half-wavelength  cavity  made 
from  a  short  section  of  standard  WR-90  X-band  waveguide.  This  brought 
about  several  problems.  We  were  attempting  to  measure  a  large  variety 
of  materials  (ceramics,  plastics,  epoxies,  laminates,  composites,  mix¬ 
tures,  etc.),  and  many  were  lossy  enough  to  lower  the  Q  below  that 
necessary  for  good  measurements.  In  other  words,  good  sharp  cavity 
resonances  on  which  to  base  the  calculations  could  not  be  obtained. 

The  high  mode  density  was  also  a  serious  problem.  The  dielectric 
constant  calculated  from  the  resonant  frequency  and  the  dimensions  of 
the  dielectric-tilled  cavity  is 


9 


NWC  TP  6643 


where  c  is  the  velocity  of  light  in  a  vacuum;  f0  is  the  resonant  fre¬ 
quency;  a,b,d  are  the  cavity  width,  height,  and  length,  respectively, 
while  m,n,p  are  integers  representing  the  mode  order  (number  of  equiva¬ 
lent  half-wavelengths  in  the  x,y,z  directions). 

Empty,  the  cavity  resonated  in  the  TE101  mode  at  resonant  fre¬ 
quency  fe.  Filled  with  dielectric,  however,  the  101  resonance  drops 
to 


✓  e' 


Since  most  radome  materials  have  e'  >  2,  filling  the  cavity  with 

dielectric  moved  the  101  mode  out  of  band  while  bringing  the  102  mode 
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within  band.  Unfortunately,  it  also  brings  others  such  as  the  TE2qj, 
TE.ij,  and  TM ^ ^  modes  within  band  (8  to  12  gigahertz). 

The  cavity  symmetry  was  such  that  normally  most  of  these  modes 
(excluding  the  TE1Q2)  should  not  be  excited.  Unfortunately,  many  of 
the  materials  being  measured  were  not  strictly  homogeneous  (ceramics, 
experimental  composites,  mixtures),  and  the  extraneous  modes  did  some¬ 
times  appear,  apparent  symmetry  notwithstanding. 

A  third  difficulty  involved  the  coupling  irises.  Equations  such 
as  Equation  2-1  are  derived  on  the  assumption  of  perfectly  conducting 
cavity  walls.  Coupling  into  the  cavity  is  usually  accounted  for  by 

assuming  a  cavity  equivalent  resonant  lumped  circuit  with  perfect 
transformers  supplying  the  coupling.  Actually,  an  iris  appears  to  sup¬ 
ply  some  shunt  susceptance  across  the  cavity,  lowering  the  resonant 
frequency.  This,  in  turn,  implies  that  the  electrical  cavity  length  is 
somewhat  greater  than  its  actual  physical  length.  Another,  more  exact, 
way  of  treating  this  problem  is  to  observe  that  a  single  mode  cannot  be 
used  to  match  exactly  the  boundary  conditions  at  the  end  walls  contain¬ 
ing  the  irises.  Instead,  the  lowest  order  mode  plus  an  infinite  sum  of 
higher  order  modes  are  required  theoretically,  although  this  sum  can  be 
terminated  at  some  point  for  all  practical  purposes.  If  the  higher 

order  modes  are  below  cutoff,  they  are  evanescent  and  decay  exponen¬ 
tially  away  from  the  iris  in  the  propagation  direction. 

Using  this  approach,  it  is  easy  to  see  that  the  presence  of  the 

dielectric  raises  the  cutoff  frequency  of  all  modes;  hence,  the  cou¬ 

pling  will  be  altered.  In  addition,  some  modes  are  no  longer  evanes¬ 
cent  and  thus  are  able  to  propagate  the  length  of  the  cavity.  That  is, 
the  coupling  coefficients  of  the  empty  and  filled  cavities  are  differ¬ 
ent;  thus,  they  have  some  effect  on  the  cavity  equivalent  electrical 
length  and  introduce  error  in  the  calculated  dielectric  constant.  This 
has  a  much  more  severe  effect  on  the  loss  tangent  measurements,  since 
the  measured  (loaded)  Q  is  a  strong  function  of  the  coupling  coeffi¬ 
cient  as  well  as  of  the  dielectric  and  cavity  wall  losses. 

A  fourth  set  of  problems  was  encountered  in  the  high-temperature 
measurements.  Measurements  in  the  TEmop  modes  require  a  close  sample 
fit,  particularly  in  the  y-  or  b-direction.  The  differential  thermal 
expansion,  particularly  between  the  plastic  and  Teflon-based  samples 
and  the  metal  cavity,  lead  to  sample  warping,  extrusion,  cracking,  and 
crumbling. 

To  deal  with  these  problems,  we  went  to  a  half-height  (b  = 
0.2  inch  instead  of  0.4  inch  in  Figure  2-2)  waveguide  X-band  cavity. 
This  decreases  substantially  tine  number  of  possibly  excited  modes  (mode 
density),  as  can  be  seen  from  Equation  2-1.  The  cavity  was  also 
lengthened  to  2.5  inches.  This  does  increase  the  overall  mode  density 
slightly,  although  not  as  severely  as  might  be  inferred  from  Equa¬ 
tion  2-1,  since  the  cutoff  frequency  of  a  rectangular  waveguide  is 
given  by  (Reference  2-1) 
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(f  ), 

c  mn 


2/e' 


(2-3) 


and  those  modes  which  are  below  cutoff  are  not  excited  in  any  case. 
However,  lengthening  the  cavity  decreases  the  iris  effects — the  change 
in  coupling  in  the  presence  of  the  dielectric.  (While  the  iris  still 
increases  the  electrical  length  of  the  cavity,  this  increase  is  a 
smaller  percentage  of  the  overall  length  leading  to  reduced  error.) 

The  long  cavity  also  allows  several  TEnop  resonances  (TE103, 
TE10t|,  and  TE105  in  the  empty  guide)  to  be  examined  across  X-band  (8  to 
12  gigahertz).  Mode  index  numbers  can  be  assigned  to  each  resonance 
and  then  varied  until  the  calculated  dielectric  constant  agrees  for  all 
of  them.  The  dielectric  ambiguity  is  thereby  removed  because  of  the 
difference  in  phase  velocity  for  each  mode  in  a  guided-wave  structure. 
This  assumes  a  dispersionless  dielectric,  which  is  usually  the  case 
across  a  single-microwave  band  for  typical  radome  materials. 

Finally,  the  cavity  was  only  partially  dielectric  filled.  The 
dielectric  occupied  a  cross  section,  as  shown  in  Figure  2-3.  The  con¬ 
figuration  of  Figure  2-3(a),  with  a  dielectric  slab  centered  lengthwise 
down  the  long  dimension  of  the  cavity,  is  used  primarily  for  accurate 
room-temperature  measurements.  The  TEmop  modes  are  employed,  and  for 
the  odd  modes  the  transcendental  equation  to  be  solved  is 


tan  d2>4mo/3)2  -  (pir/d)2  tan  d  j/e '  (  (W3)  2  -  (pir/d)2 


A  -  (3p/2fQd) 2 
it'  -  (3p/2f  d)2 


and  for  the  even  modes 


(2-4a) 
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J\  -  (3p/2f  d)2  tan  d,/e(u  /3)2  -  (p*/d)2 
o  1  o 

=  -  /e'  -  (3p/2f  d)2  •  tan  d?/(u  /3)2  -  (pir/3)2  (2-4b) 

o  o 


where  dj  =  t/2  (one-half  the  thickness  of  the  sample)  and  d2  = 
(a-t)/2.  0)o  =  2f0,  and  fQ  Is  given  in  tens  of  gigahertz  (i.e., 
10  gigahertz  means  fQ  =  1.0  in  these  units). 

The  integer  p  specifies  the  number  of  half-waves  along  the  length 
d  of  the  cavity.  The  index  m  is  determined  by  the  root  of  the  tran¬ 
scendental  equation.  The  smallest  value  of  e  in  Equation  2-4a  is  at 
m  =  1;  the  next  value  is  found  at  m  =  3,  and  so  on.  The  even  values 
come  from  the  roots  of  Equation  2-4b.  The  m  =  1  root  is  roughly  analo¬ 
gous  to  a  single  half-wavelength  across  the  width  of  the  cavity.  The 
m  =  2  root  from  Equation  2-4b  corresponds  to  a  full  wave,  etc.  The 
even  modes  are  suppressed  by  the  symmetry  of  the  situation,  while  the 
higher  order  modes  (m  =  3,  etc.)  are  below  cutoff  for  reasonable  choice 
of  thickness  t  and  dielectric  constant  in  the  range  of  typical  radome 
materials.  The  index  for  the  height  b  is  kept  to  zero  by  the  half¬ 
height.  Thus,  only  the  TEjQp  modes  are  normally  dealt  with. 

For  the  elevated  temperature  measurements,  LSM  modes  were  employed 
with  configuration  (b)  of  Figure  2-3.  This  configuration  allows  free 
sample  expansion  in  the  y-direction,  while  clearance  can  be  provided  at 
the  sample  edges  for  thermal  growth  in  the  x-  and  z-directions.  The 
electric  field  close  to  the  cavity  walls  is  small  since  the  tangential 
electric  field  must  go  to  zero  there;  hence,  the  effect  of  the  clear¬ 
ances  at  the  sample  edges  should  be  negligible.  This  was  confirmed  by 
measurements  on  a  sample  of  known  dielectric  constant,  first  as  a  tight 
fit  and  then  after  trimming  to  give  clearance. 

The  transcendental  equation  to  be  solved  for  the  LSM  mode  case  is 
(Reference  2-2) 


Ah  -  (N/a)2  -  (p/d)2  tan  tJkh  -  (N/a)2  -  (p/d)2 
=  -  e/k2  -  (N/a)2  -  (p/d)2 


•  tan  (b-t)/k2  -  (N/a)2  -  (p/d)2 


(2-5) 
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Here,  k0  =  2tt/X0,  or  in  terms  of  the  tens  of  gigahertz  units  of 
Equation  2-4,  io'0/3. 

It  should  be  pointed  out  as  a  pitfall  that  for  sufficient  sample 
thickness  combined  with  dielectric  constant  and  frequency,  the  micro- 
wave  energy  can  become  "trapped"  in  the  dielectric.  The  field  exterior 
to  the  dielectric  then  decays  exponentially  or  becomes  evanescent. 
This  takes  place  whenever 


(2-6) 


(The  thickness  enters  in  through  kQ,  since  the  greater  the  dielectric 
loading,  the  lower  the  resonant  frequency,  and  the  smaller  kD).  In 
such  cases.  Equations  2-4  and  2-5  can  be  rewritten  using  the  following 
relation: 


tan  jx  =  j  tan  h  x 


Equations  2-4a  and  2-4b  become 


tan  h  d2/(p7r/d)2  -  (  o> '  /  3  )  2  tan  d  ,  /e  '  (  u'  /  3)  2  -  (pir/d)‘ 

o  *o 


(2-7) 


(/( 3p/2f  d)2  -  1) 


~  (3p/2fQd)- 


while  Equation  2-5  becomes 


( 2-8a ) 


/(3p/2f  d)2  -  1  tan  d./e'U  /3)2  -  (pw/d)2 


=  -  /e’  -  (3p/2f  d)2  tan  h  d?/Tw  / 3)2  -  (p’r/d)2  (2-8b) 

o  o 
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/k2e'  -  (Nu/a)2  -  (pn/d)2  tan  t/k2e  -  (Nn/a)2  -  ( p  n/d ) 2 


=  E'/(Nir/a)2  +  (pw/d)2  -  k2 


tan  h(b-t)/(Nir/a)  2  +  (ptr/d)2  -  k2  (2-9) 

o 


Some  work  was  done  with  the  configuration  of  Figure  2-3(c),  par¬ 
ticularly  where  only  short  samples  were  available.  Assuming  the  sample 
is  centered  lengthwise,  then  Equations  2-4(a)  and  2—4 ( b )  or,  in  the 
case  of  trapped  waves,  Equations  2-8(a)  and  2-8(b)  can  be  used  directly 
by  interchanging  the  dimensions  "a"  and  “d."  Of  course,  the  integer  p 
will  differ  as  there  will  usually  be  only  one  or  two  spatial  variations 
along  "a"  (being  just  the  guide  width),  whereas  there  will  be  several 
along  "d."  It  follows  that  usually  not  only  the  index  p  be  1  or  1  in 
configuration  (c)  in  Figure  2-3,  but  the  argument  of  the  tangent  term 
will  be  of  higher  order  (e.g.,  between  it  and  3ir/2  for  an  odd  solu¬ 
tion).  In  configuration  (a)  of  Figure  2-3,  however,  p  will  be  some 
integer  greater  than  2  to  3,  and  the  lowcot  order  argument  (between 
zero  and  tt/2,  which  means  only  the  odd-mode  solution  need  be  used)  for 
the  tangent  will  yield  the  appropriate  dielectric  constant. 

It  is  interesting  to  note  that  the  configuration  of  Figure  2-3(c) 
can  be  analyzed  directly  by  signal  flow  graph  techniques,  whereas  this 
does  not  appear  to  be  the  case  in  Figure  2-3(a).  Work  still  is  to  be 
done  to  reconcile  the  apparently  different  expressions  obtained  for 
Figure  2-3(c)  from  signal  flow  graph  analysis  and  the  consideration  of 
tne  modes  in  a  partially  dielectric-loaded  resonator  with  perfectly 
conducting  walls.  This  comparison  could  be  of  ?mportance  in  both  error 
analysis  for  this  type  of  measurement  and  in  un  rstanding  the  effects 
of  coupling  and  direction  of  energy  propagation  on  a  cavity  resonance. 

Some  measurements  were  also  taken  in  other  configurations;  for 
example,  those  similar  to  Figures  2-3(a)  and  (c),  where  the  dielectric 
slab  was  not  centered,  or  the  use  of  two  long  samples  along  ttie  wave¬ 
guide  walls.  Such  configurations  might  have  particular  applications. 
For  example,  placing  the  sample  next  to  a  side  wall  in  the  TKm0p  mode 
ensures  minimum  electric  field  in  the  sample  and,  thus,  a  smaller  shift 
in  the  resonant  frequency.  This  could  be  useful  in  measuring  the  prop¬ 
erties  of  high  dielectric  constant  materials  where  it  is  not  desirable 
to  make  a  very  thin  sample. 

Some  measurements  were  also  made  using  the  TE j 0p  circular  wave¬ 
guide  cavity  loaded  with  a  dielectric  disk.  The  transcendental 
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equations  to  be  solved  are  very  similar  to  those  oi  the  dielectric- 
loaded  rectangular  cavity.  This  approach  needs  more  development  (for 
example,  the  residual  frequency  modulation  (FM)  from  our  sweepers  made 
it  difficult  to  obtain  clean  resonances  because  of  the  very  high  Qs  of 
TEjQp  cavities)  but  does  have  promise.  For  example,  the  electric 
field  lines  in  the  TE  j  q  circular  waveguide  mode  are  closed  loops  that 
go  to  zero  at  the  walls.  This  means  very  low  wall  losses  and  very  high 
intrinsic  Q.  It  also  means  that  a  sizeable  gap  between  the  dielectric 
sample  disk  and  the  cavity  is  of  little  consequence  since  the  electric 
field  is  negligible  there. 

Determining  the  dielectric  loss  tangent  for  the  partially 
dielectric-filled  cavity  remains  an  unfinished  problem.  The  difficulty 
arises  because  of  the  differences  in  field  configurations  between  the 
partially  filled  and  the  empty  cavity.  In  a  completely  filled  cavity, 
the  field  configuration  basically  is  unaltered  from  that  of  the  empty 
cavity;  hence,  the  Q  due  to  wall  losses  remains  unchanged,  and  one  can 
write  the  classic  expressions  such  as  (References  2-3,  2-4) 

Q“ '  =  Q^1  +  Qj1  (2-10) 


with  a  high  degree  of  approximation.  Here  Q  is  the  measured  Q  of  the 
dielectric  field  cavity,  Qe  the  Q  of  the  empty  cavity,  and  Qj  the 
dielectric  Q  or  reciprocal  loss  tangent.  There  are  approximations  in 
such  an  expression,  a  major  one  being  that  the  presence  of  the 
dielectric  does  not  appreciably  alter  the  cavity  coupling.  The 
measured  Q  is  related  to  the  sources  of  loss:  dielectric  loss,  wall 

loss,  and  losses  through  the  coupling  structures  to  the  external 
circuitry. 

The  theoretical  unloaded  Q  of  a  rectangular  cavity  with  walls  of 
known  conductivity  can  readily  be  found  for  the  lowest  order  modes  as 
(Reference  2-6) 


jrn  _  2b(a2  +  d2)3^2 _ =  pK 

4R  ad (a 2  +  d2)  +  2b(a3  +  d)3  R 

s  s 


where  n  =  n 0//e  (the  characteristic  impedance  of  free  space  divided 
by  the  square  root  of  the  dielectric  constant);  Rs  is  the  conductor 
surface  resistance;  a,b,d  are  the  cavity  dimensions;  and  K  depends  only 
on  cavity  geometry. 

! 
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Unfortunately,  this  says  the  "Q"  of  a  dielectric-filled  cavity 
changes  from  the  empty  cavity  even  with  a  lossless  dielectric,  the 
ratio  being 


Q  /eR 
e  s 


(2-12) 


(2-13) 


so  that 


R  /  f 
s  /  o 


R  \  f 
se  \  oe 


(2-14) 


It  follows  that 


JL  =  el/4 


(2-15) 


The  effect  given  by  Equation  2-15  seems  to  be  commonly  ignored  in 
the  usual  application  of  Equation  2-10,  and  we  have  to  wonder,  in  terms 
of  absolute  accuracy,  how  good  many  of  the  reported  loss  tangents  mea¬ 
sured  in  resonant  systems  really  are! 

It  seems  that  Equation  2-10  should  actually  relate  the  measured  Q 
of  a  dielectric-filled  cavity  to  the  dielectric  Qj  (reciprocal  loss 
tangent)  and  the  Q  of  a  cavity  filled  with  a  lossless  dielectric  with 
the  same  real  part  of  the  dielectric  constant,  at  least  for  low- loss 
materials.  That  is,  QP-1  in  Equation  2-10  should  be  multiplied  by 

E  !/\ 

Considering  such  pitfalls  as  illustrated  in  the  foregoing  discus¬ 
sion,  the  measurement  of  dielectric  losses  in  the  partially  dielectric- 
filled  rectangular  cavity  has  further  difficulties.  The  primary 
problem  is  that  the  field  distributions  are  different  in  the 
dielectric-loaded  and  in  the  empty  cavity;  thus,  there  seems  to  be  no 
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simple  way  to  t  L n d  the  term  in  Qc~*  (or  e^**/Qe)  from  the  real 
part  of  tlie  dielectric  constant  and  the  empty  cavity  Qe.  One  would 
need  to  measure  an  identical  geometry  of  partially  filled  cavity  with  a 
lossless  dielectric  with  the  same  real  part — an  approach  that  is  not 
•  'actical. 

An  iterative  process  might  offer  a  solution.  The  real  part  of  the 
dielectric  constant  could  he  calculated  as  if  the  material  were  loss¬ 
less.  This  could  be  used  to  find  the  TE,  l,SE,  of  LSM  field  distribu¬ 
tion.  This  is  then  integrated  over  the  cavity  walls  to  find  the  equiv¬ 
alent  Qe,  with  Qj  then  calculated  from  Equation  2-10  or  its  actual 
experimental  equivalent 


(2-16) 


where  Qc  is  the  Q  due  to  losses  to  the  external  circuit  through  the 
coupling  apertures.  (Q  would  then  be  the  actual  measured  loaded  Q. ) 
Qc  perhaps  could  be  obtained  to  sufficient  accuracy  by  comparing  the 
measured  Qe  of  the  actual  empty  cavity  with  the  theoretical  empty  Q 
as  in  Equation  2-11,  assuming  the  coupling  is  independent  of  dielectric 
loading.  Another  approach  would  be  to  conceptually  break,  the  transcen¬ 
dental  equations,  such  as  2-4,  2-5,  2-8,  or  2-9,  into  real  and  imagi¬ 
nary  parts  with  e  complex  and  attempt  to  solve  these  directly  in  terms 
of  the  experimental  measurements. 

The  complexity  of  the  problem  perhaps  can  be  underscored  in  the 
case  of  wave  trapping,  where  the  waves  in  the  empty  portions  of  the 
cavity  become  evanescent.  The  wall  Losses  must  then  become  vanishingly 
small  in  the  empty  part  of  the  cavity,  with  the  wall  losses  substan¬ 
tially  confined  to  those  parts  of  the  walls  in  direct  contact  WLth  the 
dielectric.  In  the  case  of  a  low-loss  dielectric,  the  measured  loaded 
Q  could  conceivably  be  larger  than  that  of  the  empty  cavity. 

In  summation,  the  problem  of  measuring  accurate  loss  tangents  has 
yet  to  be  satisfactorily  resolved  for  the  partially  dielectric-filled 
rectangular  cavity  and  awaits  further  efforts. 


EXPERT  MEN TA L  RESULTS 


Figure  2-4  illustrates  the  manner  in  which  the  data  are  reduced. 
The  sample  is  inserted  into  the  cavity  and  the  resonant  frequencir 
(four  in  this  example)  are  measured  across  X-band.  This  example  has 
the  sample  configuration  of  Figure  2-3(a).  The  transcendental 
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Equations  2-4  or  2-8  are  used  along  with  the  known  cavity  and  sample 
dimensions  to  calculate  the  dielectric  constant(s)  for  different 
numbers  of  half  wavelengths  m  (expressed  as  p  in  the  equations)  along 
the  cavity  length.  Figure  2-4  shows  plots  for  each  resonant  frequency 
of  the  calculated  dielectric  constant  versus  a  normalized  index,  where 
r  is  the  order  number  of  the  resonant  frequency:  r  =  1  for  f01,  r  =  2 
for  f  q  2 »  etc.  In  this  example,  the  plots  for  the  different  resonant 
frequencies  cross  for  m-r-1  =  3,  and  it  follows  that  this  is  the 
correct  unambiguous  dielectric  constant  for  this  glass  window  material 
used  in  a  high-speed  antiradiation  missile  (HARM)  application.  The 
discrepancies  (standard  deviation  0.008)  between  the  calculated  values 
are  due  to  various  experimental  errors,  while  the  curvature  in  the 
plots  is  a  result  of  waveguide  dispersion. 

This  setup  has  been  used  to  measure  a  large  variety  of  materials. 
Some  typical  room-temperature  results  are  given  in  Table  2-1. 


TABLE  2-1.  Representative  Measured  Dielectric  Constants  From 
Partially  Dielectric-Filled  Rectangular  X-Band  Cavity. 


Dielectric 

Material 

Material 

constant 

source 

Comments 

Avcoat  8029 


For  NWC,  Code  3242. 
Contains  reinforcing 
fibers;  difficult  to 
machine . 


Teflon 

Laminate 


SB-129  Propellant 


Mullite  Vat  21 


Space  shuttle  tile 


Commercial  Control. 


Naval  Air 
Development 
Center 
(NADC) 


Naval 
Research 
Laboratory 
(  NRL) 

Rockwell 

International 

Corp. 


A.  Lee,  Code  60613, 
NADC.  Very  tough 
material;  anisotropic 
average. 

Lossy;  measurement 
accuracy  poor. 

Code  6360,  NRL. 
Apparently  homogene¬ 
ous  . 

For  NWC,  Code  3522. 
Fragile;  crumbles; 
pure  Si02;  90%  air. 


Stycast  HiK  =  12 


11.70 


Emerson  and 
Cummi ngs 


Control. 
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Anisotropic  materials  also  play  a  role  in  radome  construction,  and  this 
system  can  be  used  to  measure  such  materials  in  the  TEm0p  modes 
provided  the  electric  field  is  aligned  along  a  principal  axis  of  the 
dielectric.  Some  results  of  such  measurements  are  given  in  Table  2-2. 


TABLE  2-2.  Representative  Measured  Dielectric  Constants 
for  Anisotropic  Materials  Using  Partially  Dielectric- 


Filled  Rectangular  X-Band  Cavity. 


Material 

e 

X 

e 

z 

3M  Company,  300  board 

2.83 

2.41 

Rogers  Corp.,  RT  Duroid 

D05880  (  17SK.S) 

2.  19 

2.29 

NADC  experimental  laminate 

3.07 

2.66 

Duroid  5870M 

2.34 

2.42 

This  method  can  be  adapted  to  measure  certain  types  of  inhomogene¬ 
ous  materials  by  cutting  several  small  samples  and  measuring  the  per¬ 
mittivity  of  each  piece,  using  the  configuration  of  Figure  2-3(c). 
Homogeneity  is  an  important  problem  with  radome  and  microwave  dielec¬ 
tric  materials,  which  are  mixtures,  laminates,  fiber  reinforced,  or 
ceramic.  Maintaining  even  distributions,  desired  volume  fractions,  and 
densities  poses  difficult  problems.  Table  2-3  gives  some  measured 
results  on  inhomogeneous  materials.  Omohundro  Corporation  products 
NAFA-5030  and  NAFA-2575-l  are  syntactic  foams  made  using  hollow  alumina 
microspheres,  while  the  lens  material  is  a  polyfoam  fillet!  with  alumi¬ 
num  flakes.  As  in  aluminum  spray  paint  (at  least  the  spray  paints  we 
have  checked),  the  aluminum  flakes  are  evidently  coated  with  insulating 
materials  to  prevent  f lake-to-f lake  electrical  contact.  Careful  visual 
examination  indicated  that  the  inhomogeneity  was  caused  by  an  uneven 
distribution  of  the  aluminum  particles. 

Finally,  of  course,  the  cavity  can  be  used  to  measure  dielectric 
properties  at  elevated  temperatures  up  to  800°F  in  the  LSM  modes. 
Hindsight  tells  us  that  free-space  methods  are  preferable  for  high- 
temperature  measurements  if  a  wide  variety  of  materials  is  to  be  con¬ 
sidered.  If  cavities  (or  waveguides)  are  to  be  used,  the  thermal 
expansion  of  both  the  sample  and  the  waveguide  or  cavity  must  be  con¬ 
sidered,  both  from  the  electromagnetic  and  the  structural  points  of 
view.  A  sample  that  warps  or  splits  will  probably  invalidate  the  EM 
analysis,  and  one  that  deposits  a  chemical  coating  on  a  cavity  interior 
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TABLE  2-3.  Representative  Measurements  (Relative 
Dielectric  Constant)  of  Inhomogeneous  Materials. 


Sample 

number 

Omohundro 

NAFA-5050 

Omohundro 

NAFA-2575-1 

Microwave 

lens 

material 

1 

1.95 

2.56 

2.69 

2 

1.96 

2.52 

3.25 

3 

2.00 

2.50 

2.95 

4 

1.88 

2.52 

2.88 

5 

1.89 

2.57 

3.64 

6 

1.89 

2.46 

2.85 

7 

1.94 

2.43 

2.82 

8 

1.90 

— 

3.21 

9 

— 

— 

2.82 

10 

— 

— 

3.56 

11 

— 

— 

3.48 

12 

— 

— 

3.03 

13 

— 

— 

2.69 

14 

— 

— 

2.89 

15 

— 

— 

3.46 

16 

— 

— 

2.99 

17 

— 

— 

2.82 

18 

— 

— 

2.87 

Average 

1.93+0.04 

2.51 ±0.05 

3.05+0.31 

will  probably  at  least  ruin  the  Q-value.  For  example,  Figure  2-5  shows 
some  measured  results  for  Duroid  5870M.  Here  the  dielectric  constant 
is  calculated  from  the  cavity  resonant  frequencies,  corrected  for  the 
thermal  expansion  of  the  cavity,  as  determined  by  experimental  measure¬ 
ment  of  the  resonant  frequencies  versus  temperature.  There  are  two 
curves:  (1)  dielectric  constant  calculated  without  including  the 

sample  thermal  expansion  and  (2)  decreasing  dielectric  constant  with 
increasing  temperature  calculated  with  the  sample  thermal  expansion. 

In  this  case,  the  correction  for  the  expansion  of  the  beryllium- 
copper  cavity  is  relatively  insignificant  compared  to  that  brought 
about  by  the  large  thermal  expansion  coefficient  of  Duroid.  The  ther¬ 
mal  expansion  data  used  are  shown  in  Figure  2-6.  In  this  case,  the 
reinforcing  fibers  lie  in  the  plane  of  the  sample,  and  the  thermal 
expansion  coefficient  normal  to  the  plane  was  used.  In  the  case  of 
materials  with  such  a  large  thermal  expansion,  the  calculated  dielec¬ 
tric  constant  can  be  no  better  than  the  thermal  expansion  coefficient 
used.  Figure  2-7  shows  some  data  obtained  from  the  Rogers  Corporation 
for  both  5650M  and  5870M  Duroids.  While  the  curves  of  Figure  2-7(a) 
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FIGURE  2-6.  Thermal  Expansion  for  Duroid  5870M 
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Since  the  radome  expands  with  increasing  temperature,  one  could 
argue  that  the  "effective"  dielectric  constant,  which  would  be  repre¬ 
sented  by  the  curve  of  increasing  dielectric  constant  (with  tempera¬ 
ture)  in  Figure  2-5,  is  the  quantity  of  most  interest.  Using  this 
value,  perhaps  the  radome  could  be  treated  as  having  a  constant  wall 
thickness  but  changing  effective  dielectric  constant  to  a  reasonable 
degree  of  approximation.  In  an  actual  missile  flight,  the  radome  will 
not  come  to  some  uniform  temperature  but  rather  will  sustain  tempera¬ 
ture  gradients  that  can  be  very  substantial.  On  a  low  thermal  conduc¬ 
tivity  material  such  as  5870M  during  a  short  missile  flight,  only  a 
fairly  thin  outer  layer  may  undergo  a  large  change  in  temperature.  The 
use  of  an  "effective"  dielectric  constant  while  holding  the  radome 
dimensions  constant  with  an  "N-layer"  analysis  seems  attractive, 
although  its  adequacy  remains  to  be  shown. 


NWC  TP  6643 


REFERENCES 

2-1.  R.  E.  Collin.  Foundations  for  Microwave  Engineering.  New  York., 
McGraw-Hill,  1966.  P.  323. 


2-2. - .  Field  Theory  of  Guided  Waves.  New  York,  McGraw-Hill, 

1966.  Pp.  224-32. 

2-3.  - .  Foundations  for  Microwave  Engineering.  New  York, 


McGraw-Hill,  1966.  P.  315. 

2-4.  "High  Frequency  Dielectric  Measurement,"  in  Proceedings  of  a 
Tutorial  Conference  on  Measurement  of  High  Frequency  Dielectric 
Properties  of  Materials,  National  Physical  Laboratory,  Tedding 
ton,  U.K.  27-29  March  1972.  J.  Chamberlain  and  G.  W.  Chantry 
eds.  Place,  PIC  Science  and  Technology  Press,  Ltd.,  1973.  Pp 
12-13. 

2-5.  S.  Ramo  and  J.  R.  Whinnery.  Fields  and  Waves  in  Modern  Radio 
New  York,  John  Wiley  and  Sons,  1953.  Pp.  424-26. 


NWC  TP  6643 


MEASURING  COMPLEX  PERMEABILITY  AND  DIELECTRIC  CONSTANT  OF 
ISOTROPIC  MATERIALS  IN  A  COAXIAL  TRANSMISSION  LINE 


A  Hewlett-Packard  (HP)  8410B  automatic  network  analyzer  is  used  to 
measure  the  complex  reflection  Sll  and  transmission  S21  coefficients  of 
a  homogeneous  sample  that  is  placed  in  a  7-millimeter  air  line  as  shown 
in  Figure  3-1.  The  air  line  connectors  and  air  line  are  assumed  to  be 
lossless.  The  reflection  coefficients  Sll  and  S22  of  the  air  line  con¬ 
nectors  are  assumed  to  be  zero.  The  line  length  from  the  measurement 
plane  to  the  sample  is  the  only  parameter  that  is  taken  into  account. 
This  system  allows  measurements  over  a  frequency  range  of  1  to 
12.4  gigahertz.  The  frequency  range  is  limited  by  the  HP-8746  S- 
parameter  test  set  being  used  in  the  current  setup. 


OUTER  CONDUCTOR 


CENTER  CONDUCTOR 
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The  HP-8410B  automated  network  analyzer  system  provides  removal  of 
the  systematic  errors  caused  by  mismatch  and  leakage  in  the  test  setup, 
isolation  between  the  reference  and  test  channels,  and  the  frequency 
tracking  errors  between  test  channels.  These  errors  are  removed  by  use 
of  an  8-term  error  model  (Reference  3-1).  With  this  accuracy- 
enhancement  software,  the  amplitude  and  phase  of  the  scattering  parame¬ 
ters  can  be  measured  to  within  0.1  decibel  and  1  degree,  respectively. 
The  sample  holder  is  a  7-millimeter  coaxial  air  line  with  k PC-7  con¬ 
nectors.  The  sample  is  accurately  positioned  in  the  air  line  by  means 
of  a  probe  with  a  micrometer  dial  attached.  Sllabc  and  S21abc  of  the 
sample  and  sample  holder,  respectively,  are  measured,  and  Sllb  and  S12b 
of  the  sample  are  given  by  Equations  3-1  and  3-2  below: 


Sllb  = 


S 1 labc 
S2 la  S 1 2a 


(3-1) 


Sllb  is  the  reflection  coefficient  of  the  sample.  Sllabc  is  the  mea¬ 
sured  rejection  coefficient  of  the  sample,  air  line,  and  its  connec¬ 
tors  (see  Figure  3-2).  S21a  and  SI  2a  are  the  forward  and  reverse 


AIR  LINE  AND 
CONNECTOR 


SAMPLE 


AIR  LINE  AND 
CONNECTOR 
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S12b 


S21abc 
S12a  S12c 


(3-2) 


S12b  and  S21abc  are  the  reverse  transmission  and  the  measured  forward 
transmission  coefficients,  respectively,  of  the  sample.  S12a  and  S12c 
are  the  reverse  transmission  coefficient  of  the  connectors  a  and  c  and 
the  air  line  up  to  the  sample  interface. 


S12c 


lc  is  the  electrical  line  length  from  the  sample  interface  to  the 
Port  2  measurement  plane. 

Since  we  now  have  the  complex  reflection  coefficient  SI  lb  and  the 
complex  transmission  coefficient  S12b  of  the  sample,  we  can  now  compute 
the  following  complex  equations  that  will  give  us  the  complex  dielec¬ 
tric  constant  cc  and  complex  permeability  uc*  The  derivations  of 
the  following  equations  are  presented  in  Reference  3-2. 

\  j  =  SI  lb  +  S12b  (3-3) 


V2  =  S 1 2b  -  Sllb 


(3-4) 


(1  ~  VjV2) 
Vl  -  v2 


(3-5) 


r  =  X  ±  A2  -  1 


(3-6) 


Here  r  is  the  reflection  at  the  boundary  of  the  sample.  Root  r  is 
chosen  so  that  Irl  <  1. 
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(1  +  r)‘ 


C  H2 

“  ~2  d 

to 


L  i  +  .  i  \2 

(  in  -  +  j  arg  - \ 

l  M  lzl j 


(3-8) 


(3-9) 


Here  w  =  2nv  and  d  is  Che  thickness  of  the  sample.  The  arg  1/Z  = 
2im  cos-*  (Z'//Z|),  Z'  is  the  real  part  of  Z,  and  n  is  an  integer  that 
defines  an  ambiguity  that  occurs  every  2ir  radians.  From  Equations  3-8 
and  3-9,  the  complex  dielectric  constant  is  given  by 


(3-10) 


/C  jC  2 


(3-11) 


ERROR  ANALYSIS 


After  measuring  a  few  samples,  including  air,  it  became  apparent 
that  variations  in  the  dielectric  constant  were  systematic  with  fre¬ 
quency.  These  systematic  errors  were  caused  by  a  load  mismatch  and 
were  not  taken  into  account  when  calibrating  the  network  analyzer  (the 
8-term  error  model  does  not  include  isolation  and  load-match  errors). 
The  12-term  error  model  requires  that  all  four  S-parameters  be  mea¬ 
sured,  but  it  provides  for  the  best  measurement  accuracy  of  magnitude 
and  phase.  In  the  near  future,  the  12-term  error  model  will  be  incor¬ 
porated  to  determine  if  the  dielectric  constant  accuracy  can  be 
improved  over  the  technique  that  is  presently  being  used. 

To  reduce  measurement  time,  the  present  technique  is  to  place  a 
10-decibel  pad  at  the  Port  2  measurement  plane  to  reduce  the  voltage 
standing  wave  ratio  (VSWR).  The  measured  scattering  parameters  of  an 
empty  air  line  are  then  vectorlally  subtracted  from  the  measurements 
performed  with  the  sample  in  the  air  line.  When  measuring  air  with 
this  technique,  a  dielectric  constant  of  1.000  is  calculated  with  a 
standard  deviation  of  0.008.  For  Rexolite,  the  dielectric  constant  is 
2.544  with  a  standard  deviation  of  0.047.  Figures  3-3  and  3-4  are 
plots  of  these  measurements. 


XS 
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EMPTY  AIR  LINE 

AVERAGE  DIELECTRIC  =  1  000 

STANDARD  DEVIATION  =  0  005 


FIGURE  3-3.  Measured  Dielectric  Constant  of  Air  in  the  Empty  Air 
Line. 


To  determine  the  effect  of  measurement  errors  SI  1  and  S21  on  the 
dielectric  constant  and  permeability,  known  errors  are  imposed  on  cal¬ 
culated  values  of  Sll  and  S12  and  new  values  of  dielectric  constant  and 
permeability  are  calculated.  At  this  time,  the  errors  study  is  just 
beginning.  If  the  measurement  system  is  the  major  source  of  error, 
then  additional  improvements  can  be  made,  such  as  employing  a  phase- 
locked  radio  frequency  (RF)  source  to  prevent  harmonic  skip.  In  this 
case,  a  different  harmonic  is  selected  by  the  network  analyzer  for  the 
same  test  frequency  between  the  time  of  calibration  and  measurement. 
Sample  position  error  should  be  held  to  less  than  the  1- degree  phase 
measurement  error,  which  would  be  2.6  mils  at  12.6  gigahertz. 
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REXOLITE  CENTERED  IN  AIR  LINE.  0  241  INCH 

AVERAGE  DIELECTRIC  *  2  544 
STANDARD  DEVIATION  =  0.047 


MHz/DIV 

FIGURE  3-4.  Measured  Dielectric  Constant  of  Rexolite  in  Coaxial  Air 
Line . 


Waveguide  modes  can  propagate  in  high-dielectric  or  permeability 
materials.  The  cutoff  frequency  for  the  TE j j  mode  is  given  by 


V 

f  =  0.293  - 
c  a 


where  V  is  the  velocity  of  light  through  the  medium  given  by 
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C  is  the  speed  of  light,  and  er  is  the  relative  dielectric  constant 
of  the  medium.  Thus,  for  air  the  frequency  at  which  a  TEjj  mode  can 
propagate  in  a  7-millimeter  air  line  is  approximately  12.5  gigahertz. 
The  waveguide  mode  can  exist  only  in  the  sample  itself,  since  the  coax¬ 
ial  line  will  support  only  a  transverse  electromagnetic  (TEM)  wave  at 
this  frequency.  Since  the  waveguide  modes  cannot  be  propagated  in  the 
air  line,  the  energy  is  dissipated  in  the  sample  and  makes  the  sample 
appear  lossier  than  it  really  is.  The  strength  of  the  waveguide  mode 
is  dependent  upon  the  efficiency  of  the  conversion  of  the  TEM  wave  to  a 
waveguide  mode.  This  efficiency  can  be  minimized  by  ensuring  that  the 
sample  faces  are  perpendicular  to  the  longitudinal  axis  of  the  sample 
and  that  the  faces  are  smooth.  Another  source  of  error  that  needs  to 
be  investigated  is  the  effect  of  a  small  air  gap  between  the  sample  and 
the  inner  and  outer  conductors  of  the  air  line. 
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FREE-SPACE  BRIDGE  DIELECTOMETER 

A  free-space  bridge  is  well  suited  for  measuring  the  dielectric 
properties  of  materials  at  high  temperature,  since,  ideally,  the  sample 
can  be  heated  without  producing  effects  caused  by  heating  the  measure¬ 
ment  apparatus.  In  fact,  the  sample  can  be  heated  beyond  the  damage 
threshold,  allowing  destructive  testing.  The  free-space  bridge  offers 
other  advantages  over  cavity-type  dielectric  measurements.  For 
instance,  the  free-space  bridge  is  generally  used  with  simple  flat 
panel  samples;  whereas,  for  accurate  results,  closed-cavity  measure¬ 
ments  require  machining  the  sample  to  very  close  tolerance.  Also,  the 
polarization  and  incidence  angle  of  the  impinging  wave  are  widely  vari¬ 
able  in  the  free-space  bridge.  This  is  especially  useful  for  studying 
anisotropic  materials. 

The  free-space  bridge  does  pose  some  specific  disadvantages  com¬ 
pared  with  a  cavity  dielectometer.  While  very  precise  measurements  can 
be  obtained  in  the  free-space  bridge  dielectometer ,  its  inherent  accur¬ 
acy  is  degraded  by  the  inevitable  presence  of  undesired  reflections. 
The  free-space  measurement  can  track  very  small  changes  in  the  dielec¬ 
tric  constant  with  change  in  temperature  and  can  detect  small  differ¬ 
ences  between  similar  samples.  However,  for  absolute  accuracy,  these 
difference  measurements  must  be  referenced  to  a  measurement  with  better 
inherent  accuracy,  usually  a  resonant  cavity-tvpe  measurement. 

At  NWC,  a  free-space  bridge  dielectometer  has  been  constructed  to 
measure  the  complex  transmission  coefficient  of  flat- panel  samples  at 
35  gigahertz.  The  complex  permittivity  of  the  sample  is  calculated 
from  the  transmission  measurements  by  iteration  using  a  digital  com¬ 
puter.  An  array  of  high-intensity  heat  lamps  is  used  to  heat  the 
sample;  bulk  temperatures  up  to  1000°C  have  been  obtained  thus  far. 

A  block  diagram  of  the  apparatus  is  given  in  Figure  4-1.  The 
microwave  beam  is  directed  through  the  sample  by  a  pair  of  lens-focused 
horns.  For  the  current  setup,  the  lens  is  focused  to  a  3-dec i be  1 -spot 
size  of  1.5  wavelengths  at  an  18-inch  focal  distance.  Because  of  spot- 
focusing  the  beam,  samples  only  a  few  inches  in  diameter  can  be  used 
without  significant  errors  caused  by  diffraction  at  edges  of  the  sam¬ 
ple.  The  receive  horn  sits  upon  a  carriage  which  moves  the  horn  along 
and  across  the  beam;  micrometer  screws  read  the  position  of  the  horn  to 
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FLGURE  4-1.  Block  Diagram  of  the  Free-Space  Bridge 
Die lectometer. 


0.G001  inch.  The  reference  arm  of  the  bridge  contains  a  phase  shifter 
and  a  variable  attenuator  that  are  used  to  precisely  "null"  and  balance 
the  bridge  prior  to  inserting  the  test  sample. 

The  signal  from  the  free-space  arm  of  the  bridge  is  combined  with 
the  signal  from  the  reference  arm  in  a  tuned  hybrid.  The  combined  sig¬ 
nals  reach  a  crystal  detector  producing  a  direct  current  voltage  that 
depends  on  the  relative  phase  and  amplitude  of  the  signals  from  the  two 
bridge  arms.  A  lock-in  amplifier  measures  the  detected  signal  voltage 
and  provides  greatly  increased  sensitivity  over  the  use  of  a  direct 
current  voltmeter,  permitting  the  measurement  of  a  very  small  power 
imbalance  on  the  bridge.  This  enhances  the  precision  of  the  measure¬ 
ment  and,  in  principle,  allows  the  measurement  of  transmission  loss 
through  very  low-loss  materials. 

The  RF  signal  source  is  a  10-milliwatt  klystron,  operating  at 
35  gigahertz,  that  is  phase-locked  to  a  stable  5-megahertz  quartz  crys¬ 
tal  oscillator.  Without  phase- locking  the  source,  precision  bridge 
measurements  would  be  limited  by  frequency  drift.  The  35-gigahertz 
signal  is  square-wave  modulated  to  provide  a  signal  for  the  lock-in 
amplifier. 

The  sample  heater  consists  of  a  set  of  ten  high-intensity  incan¬ 
descent  lamps  with  collimating  parabolic  reflectors.  The  lamps  are 
arrayed  in  a  circular  arc  on  either  side  of  the  sample  to  focus  their 
!  output  onto  a  narrow  strip  on  the  sample  surface.  At  full  power,  each 
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lamp  produces  an  output  heat  flux  of  approximately  1  kilowatt  on  a  1.5— 
by  6-inch  aperture  with  peak  output  at  ~l-micron  wavelength.  The 
lamps  are  cooled  by  a  recirculating  water  cooler. 

The  sample  temperature  is  measured  by  a  thermocouple  placed  below 
the  surface  of  the  sample.  A  programmable  temperature  controller, 
which  is  calibrated  for  the  thermocouple,  provides  control  of  the  bulk 
temperature  of  the  sample.  The  controller  senses  the  temperature  of 
the  sample  and  sends  a  signal  voltage  to  a  silicon-controlled  rectifier 
(SCR)  power  controller  that  adjusts  the  lamp  power  to  maintain  a  preset 
sample  temperature  or  to  change  the  temperature  on  a  programmed  time 
schedule. 

Figure  4-2  is  a  photograph  of  the  NWC  free-space  bridge  dielectom- 
eter  taken  with  the  heat  lamps  turned  on,  producing  the  diffuse  bright 
spot  in  the  center  of  the  picture.  The  transmit  and  receive  horns  are 
mounted  on  separate  supports  assembled  from  concrete  blocks.  The 
sample  holder  and  heat  lamp  arrays  are  mounted  together  on  a  table 
that  also  supports  the  temperature/power  controllers  and  a  vent  hood. 
The  purpose  of  the  vent  hood  is  to  draw  away  some  of  the  waste  heat 
produced  by  the  lamps  and  to  remove  any  gases  that  might  be  evolved  by 
the  hot  sample.  Also  shown  in  Figure  4-2  is  the  rack  that  carries  the 
electronic  instrumentation  for  the  dielectometer . 


FIGURE  4-2.  NWC  Free-Space  Bridge  Dielectometer. 
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Initial  calibration  of  the  dielectometer  is  made  by  setting  the 
phase  shifter  and  attenuator  in  the  bridge  reference  arm  so  that  the 
bridge  is  precisely  nulled  and  balanced,  i.e.,  the  output  of  the  detec¬ 
tor  is  zero  without  the  sample  in  place.  With  the  sample  in  place,  the 
bridge  is  renulled  at  each  temperature  by  adjusting  the  separation  of 
the  focused  horns  to  obtain  a  minimum  of  the  detected  voltage.  The 
phase  shift  through  the  sample  is  measured  from  the  displacement  of  the 
micrometer,  which  moves  the  receive  horn  along  the  beam  axis.  This 
method  of  phase-shift  measurement  offers  the  best  possible  precision; 
it  gives  an  accurately  symmetrical  interference  minimum  and  restores 
the  electrical  distance  between  the  horns  to  the  same  value  as  the  ini¬ 
tial  calibration  (plus  some  integer  number  of  hal f-wavelengths ) . 

Transmission  loss  through  the  sample  is  determined  by  measuring 
the  power  imbalance  between  the  two  bridge  arms  when  the  bridge  has 
been  adjusted  to  give  a  detected  interference  minimum  with  the  sample 
in  place.  The  power  imbalance  can  be  read  directly  from  the  variable 
attenuator  in  decibels  of  attenuation  required  in  the  reference  arm  to 
zero  the  detector  output. 

The  magnitude  of  the  sample  voltage  transmission  coefficient  can 
be  determined  from  detector  voltage  readings  without  the  need  of  detec¬ 
tor  calibration,  provided  that  the  detector  exhibits  proper  square-law 
behavior  over  the  dynamic  range  of  the  measurement.  Using  superposi¬ 
tion  of  signals  from  the  bridge  reference  and  test  arms  for  interfer¬ 
ence  maximum  and  interference  minimum,  the  magnitude  t  of  the  sample 
voltage  transmission  coefficient  can  be  obtained  in  terms  of  detector 
output  voltage  levels.  The  quantity  t  is  given  by  the  following 
expression: 


t  = 

1  +  /P 


(4-1) 


where  P  is  ^max^min  —  the  ratio  of  minimum  to  maximum  detected 
power  as  the  relative  bridge  phase  is  rotated  360  degrees  with  the  sam¬ 
ple  in  place.  For  a  good  square-law  detector,  P  is  the  ratio  of  the 
maximum  and  minimum  detector  voltages. 

Except  in  the  case  of  normal  incidence  or  for  very  thin  samples, 
displacement  of  the  microwave  beam  by  refraction  through  the  sample  is 
important  to  the  measurement  of  the  transmission  loss  through  the  sam¬ 
ple.  For  this  reason,  the  receive  horn  can  be  moved  laterally  to  the 
beam,  using  a  micrometer  screw  to  offset  the  deflection  of  the  beam  by 


NWC  TP  6643 


the  sample.  For  a  sample  of  thickness  d  and  relative  dielectric  con¬ 
stant  er,  the  beam  deflection  as  a  function  of  incidence  angle  A  is 
given  by 


s  =  d  sin  0  1  - 


cos  0 


'er  -  sin  0 


(4-2) 


Since  the  determination  of  the  relative  dielectric  constant  er  is 
nearly  independent  of  the  sample  loss,  it  can  be  approximately  measured 
without  regard  to  beam  deflection  and  then  employed  in  the  equation 
above  for  making  a  refined  measurement. 

Having  measured  the  phase  shift  AL  through  the  sample  and  the 
voltage  transmission  magnitude  t,  the  relative  dielectric  constant  and 
loss  tangent  of  the  test  material  are  approximately  given  by  the 
following  expressions  (Reference  4-1): 


e  =  1  +  2 
r 


cos  0  + 


(4-3) 


tan  6 


X£n(t)/er  -  sin2  0 


(4-4) 


Here  X  is  the  free-space  wavelength  and  0  is  the  plane-wave  angle  of 
incidence  on  the  sample.  These  equations  are  valid  for  either  linear 
polarization,  but  they  are  only  approximate  because  they  do  not  account 
for  the  effect  of  internal  reflections  at  the  sample  surfaces.  For  the 
case  in  which  the  electric  field  polarization  is  parallel  to  the  plane 
of  incidence  (parallel  polarization),  the  error  produced  by  use  of 
these  equations  is  generally  smaller  than  for  perpendicular  polariza¬ 
tion.  In  fact,  for  parallel  polarization,  the  equations  are  exact  at 
an  angle  of  incidence  equal  to  the  Brewster  angle, 


0g  =  tan" 


(4-5) 


in  which  case  the  interface  reflection  coefficient  r^  is  equal  to 
zero.  For  either  polarization,  the  maximum  measurement  error  of  AL  is 
equal  to  (Reference  4-2) 
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jL  sin  1  (r  ) 2  (4-6) 

"I  *■ 


For  arbitrary  polarization  and  angle  of  incidence,  the  interface 
ret  lection  from  the  sample  cannot  be  ignored  and  the  generalized  flat- 
panel  transmission  coefficient  must  be  employed.  The  complex  transmis¬ 
sion  coefficient  for  a  plane  wave  incident  on  a  flat  sample  in  air  can 
be  obtained  by  satisfying  boundary  conditions  and  using  Snell’s  law  at 
both  sample  surfaces  for  incident,  reflected,  and  transmitted  waves 
(Reference  4-3).  The  transmission  coefficient  thus  obtained  for  non¬ 
magnetic  materials  (p  =  p0)  is  given  by 


(1  -  r  2)exp  -j  6  d(cos  6  +  /e  -  sin2  9) 

T  = - - - 5 - - - - - -  (4-7) 

1  -  ^ 2exp(-2j 8Qd/er  -  sin2  9) 


where  d  is  the  sheet  thickness  and  9  is  the  angle  of  incidence;  B0  is 
the  free-space  wave  number  equal  to  2ir/A0.  The  relative  dielectric 
constant  of  the  sample  sheet  is  er*  For  lossy  material,  this  is 
replaced  by  the  complex  relative  dielectric  constant,  er(l  -  j  tan  6), 
in  which  tan  6  is  the  so-called  loss  tangent  or  dissipation  factor  of 
the  material.  The  interface  reflection  coefficient  or  Fresnel  coeffi¬ 
cient  r^  depends  on  the  polarization  of  the  incident,  wave.  For  elec¬ 
tric  field  parallel  to  the  plane  of  incidence  requiring  tangential  E 
and  H  fields  to  be  continuous  across  the  interface  gives  the  following 
formula  for  interface  reflection  coefficient  (Reference  4-4): 


ercos  9  -  /er  -  sin2  9 

r  =  _ E -  -r  (4-8) 

ercos  9  +  /er  -  sin2  9 

For  electric  field  polarized  perpendicular  to  the  plane  of  incidence, 
the  same  boundary  conditions  will  stipulate  the  following  relation  for 
ri  = 


cos  9  -  /cr  -  sin2  9 

r  =  -  ■■  -  (4-9) 

cos  9  +  /cr  -  sin2  9 
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Once  again,  for  the  validity  of  these  equations,  we  have  assumed  that 
the  material  is  nonmagnetic  and  surrounded  by  air  or  vacuum. 

By  the  measurement  procedures  previously  described,  we  obtain  the 
magnitude  and  phase  of  the  complex  transmission  coefficient,  given  by 
Equation  4-7,  for  the  sample  being  tested.  With  all  of  the  other 
measurement  parameters  known,  the  equation  can,  in  principle,  be  solved 
for  the  complex  dielectric  constant  of  the  sample.  Lacking  a  formula 
that  expresses  the  dielectric  constant  in  terms  of  the  measurement 
parameters,  Equation  4-7  must  be  solved  by  numerical  techniques  using  a 
digital  computer. 

The  numerical  approach  proceeds  as  follows.  Using  an  initial 
estimate  for  the  dielectric  constant,  a  test  value  for  the  complex 
transmission  coefficient  is  calculated  from  Equation  4-7.  The  test 
value  is  subtracted  from  the  measured  value  of  the  transmission  coeffi¬ 
cient  to  define  a  test  function  of  the  dielectric  constant.  The  root 
of  the  test  function,  which  we  obtain  using  an  algorithm  based  on 

Muller's  method  (Reference  4-5),  is  the  desired  complex  dielectric  con¬ 
stant  of  the  material.  The  computer  program,  which  is  run  on  an  HP- 
1000  computer,  appears  in  Appendix  A. 

Actually,  there  is  an  infinite  number  of  dielectric  constant 

values  which  will  solve  Equation  4-7  for  a  given  value  of  T.  This  does 
not  present  a  problem  if  the  initial  estimate  is  reasonably  close  to 

the  actual  value.  For  a  totally  unknown  sample,  a  very  good  initial 
estimate  can  be  obtained  by  finding  the  angle  of  incidence  that  gives 
maximum  transmission  through  the  sample  for  parallel  polarization. 
This  being  the  Brewster  angle,  the  relative  dielectric  constant  is 

obtained  from  Equation  4-5,  although  the  precision  is  limited  by  the 
precision  of  the  angle  measurement. 

Dielectric  constant  data  have  been  obtained  for  a  few  sample 
materials  at  room  temperature,  and  the  results  are  presented  in 
Table  4-1.  These  are  primarily  electromagnetic  window  materials  that 
have  been  provided  by  industry.  The  data  were  taken  over  a  range  of 
incidence  angles,  and  the  relative  dielectric  constant  given  in  the 
table  represents  an  average  of  the  data  for  each  sample,  assuming  the 
materials  are  basically  isotropic.  The  uncertainty  that  is  given  is 
equal  to  one  standard  deviation  for  the  data  set.  Table  4-2  gives 
reference  data  provided  by  the  manufacturer  or  reported  by  independent 
investigators . 
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TABLE  4-1. 

Room  Temperature 

Manufacturer 

Material 

er 

Rogers  Corp. 

Duroid  5880 

— 

- 

Rogers  Corp. 

Duroid  5870 

2.41 

+ 

0.02 

Omohundro  Corp. 

1-Cerez 

1.98 

+ 

0.02 

Omohundro  Corp. 

1-Cerez 

1.93 

+ 

0.03 

Corning  Glass  Works 

Pyroceram  9606 

5.64 

± 

0.04 

Dynasil  Corp. 

Dynasil  4000 
(fused  quartz) 

3.77 

± 

0.05 

TABLE  4 

-2.  Dielectric 

Reference  Data. 

Source 

Material 

e 

r 

tan  6 

Manufacturer 

Duroid  58803 

2.20  ±  0.02 

0.0009 

Manufacturer 

Duroid  5870a 

2.33  ±  0.02 

0.0012 

Manufacturer 

1-Cerez 

1.95  -  2.15 

0.017 

Manufacturer 

2-Cerez 

1.65 

0.003 

Reference  4-5 

Pyroceram  9606 

5.57 

0.0003 

Reference  4-6 

Dynasil  4000 
(fused  quartz) 

3.82 

0.0003 

aData  taken 

at  10  gigahertz. 

Loss  tangent  measurements  for  these  low-loss  materials  were  found 
to  vary  widely  with  incidence  angle  because  of  systematic  effects; 
thus,  the  results  are  not  reported  here.  While  the  accuracy  was  poor, 
the  sensitivity  obtained  for  loss  tangent  measurements  using  this 
apparatus  is  quite  good.  Changes  in  transmitted  power  through  the 
sample  can  easily  be  resolved  to  0.01  decibel,  making  the  apparatus 
very  useful  for  measuring  changes  in  the  value  of  the  sample  loss 
tangent  as  it  is  heated.  For  absolute  accuracy,  measured  changes  can 
be  referenced  to  baseline  data  obtained  at  room  temperature  with  a  more 
accurate  cavity-type  measurement,  such  as  the  open  resonator 
measurements  described  in  the  following  section. 

High-temperature  measurements  of  the  relative  dielectric  constant 
have  been  obtained  for  fused  quartz  up  to  a  temperature  of  900°F. 
Higher  temperatures  could  be  obtained  by  blackening  the  surface  of  the 
sample  since  quartz  is  nearly  transparent  at  the  peak  output  wavelength 
of  the  heat  lamps.  Figure  4-3  shows  a  plot  of  the  measured  results 
along  with  measurements  made  by  W.  Ho  of  Rockwell  Science  Center 
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(Reference  4-7). 
and  our  results  < 


TEMPERATURE.  °F 


FIGURE  4-3.  High-Temperature  Measurements 
of  Fused  Quartz  Sample. 


To  obtain  the  best  possible  sensitivity  and  accuracy  with  the 
apparatus,  a  number  of  possible  sources  of  error  must  be  considered: 

Mechanical  backlash  in  the  recelve-horn  positioner 

Instrument  drift  over  time 

Sample  thickness  variability 

Sidelobes  and  diffraction  at  sample  edges 

Thermal  effects  due  to  indirect  heating  of  the  apparatus 

Presence  of  standing  waves 

The  mechanical  backlash  problem  occurs  when  the  direction  is 
reversed  while  adjusting  the  positioner  to  null  the  bridge.  The  error 
can  be  as  large  as  0.001  inch  in  air  but  can  be  eliminated  almost 
entirely  by  making  the  final  adjustment  in  the  forward  direction. 
After  an  initial  2-hour  warm-up,  instrument  drift  is  insignificant  if 
the  measurement  is  carried  out  within  10  minutes  or  so  of  initial  cali¬ 
bration.  For  best  accuracy,  the  sample  should  have  uniform  thickness 
to  at  least  the  desired  precision  of  the  measurement.  However,  since 
the  measurement  is  actually  an  average  over  the  surface  of  a  spot 
illuminated  by  the  Gaussian  beam,  an  appropriate  average  value  taken 
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for  the  sample  thickness  should  give  satisfactory  results.  While 

antenna  side  lobes  are  certainly  present,  they  produce  no  measurable 
effect  in  the  apparatus  unless  they  are  scattered  by  obstacles  placed 
in  certain  "hot  spots"  near  the  antennas.  We  have  found  that  diffrac¬ 
tion  at  the  edges  of  the  sample  is  insignificant  if  the  sample  presents 
a  cross  section  of  as  little  as  3  inches  across  the  microwave  beam. 

Indirect  heating  of  the  antenna  focusing  lenses  can  produce  large 
errors  in  a  short  time  if  thermal  isolation  is  not  provided  while  the 
sample  is  being  heated.  This  problem  is  solved  by  placing  a  shield  of 
refractory  material  in  front  of  each  lens  while  the  heat  lamps  are  on 
and  briefly  turning  off  the  lamps  while  a  measurement  is  being  made. 

Although  there  are  isolators  ahead  of  the  detector  and  in  each  arm 

of  the  bridge,  there  is  a  standing-wave  ratio  on  the  order  of  1.2:1 
between  the  two  focused  horns.  Additional  reflections  anywhere  on  the 
bridge  will  produce  error  dependent  on  the  electrical  separation  of  the 
reflections.  At  this  point,  it  is  impractical  to  try  to  further  reduce 
the  reflections  in  the  system,  but  it  may  be  possible  to  quantify  the 
effects  of  the  standing  waves,  using  an  appropriate  model  for  the 
system,  and  subsequent  ly  subtract  these  effects  from  the  measurement. 
A  model  has  been  developed  that  may  prove  useful  for  enhancing  the 
accuracy  of  the  free-space  die lectometer .  A  graphical  exposition  of 
the  model  is  presented  in  Appendix  B. 

At  the  present  time,  the  most  basic  limitation  to  the  accuracy  of 
measurements  in  the  lens-focused,  free-space  bridge  is  the  result  of 
the  Gaussian  character  of  the  wave  front  that  is  incident  on  the 
sample.  The  dielectric  constant  is  calculated  from  the  measured 
transmission  coefficient  according  to  Equation  4-7,  which  is  actually 
derived  for  the  case  of  an  incident  plane  wave.  Since  the  phase  and 
amplitude  are  complicated  functions  of  position  in  the  incident  beam. 
Equation  4-7  gives  only  an  approximate  description  of  transmission 
through  the  sample,  and  the  errors  produced  increase  with  higher 
dielectric  constant.  At  this  point,  it  is  unclear  whether  it  may  be 
practical  or  possible  to  formulate  a  model  that  accurately  accounts  for 
the  Gaussian  character  of  the  beam  in  this  type  of  measurement.  If  it 
can  be  done,  this  would  be  a  significant  accomplishment  because  of  the 
inherent  precision,  ease,  and  flexibility  of  the  free-space  bridge 
measurement . 

Other  work  to  be  undertaken  in  the  near  future  includes  the  addi¬ 
tion  of  a  reflection-measurement  capability  to  this  apparatus.  By  mea¬ 
suring  complex  reflection  and  transmission  coefficients  of  a  sample, 
the  permeability  and  permittivity  of  the  sample  can  both  be  obtained. 
The  reflection  measurement  capability  would  also  be  useful  for  the 
study  of  microwave  and  millimeter-wave  absorbing  materials. 

Since  the  lens-focused  horn  antennas  used  in  this  apparatus  are 
essentially  frequency  independent,  the  apparatus  could  be  used  in  other 
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frequency  ranges  by  the  addition  of  the  appropriate  waveguide 
"plumbing."  Current  plans  include  testing  the  apparatus  in  Ku  band. 

The  apparatus  could  be  partly  automated  using  a  readily  available, 
high-precision  linear  positioner  under  computer  control.  By  providing 
a  digital  output  for  the  detector  voltage  to  the  computer,  the  data 
could  be  acquired  automatically  by  the  use  of  a  null-seeking  algorithm 
to  position  the  receive  horn.  This  approach  could  greatly  increase  the 
accuracy  of  the  data  by  eliminating  mechanical  backlash  error  and  by 
minimizing  sample  temperature  drift  in  a  much  faster  mt  surement  cycle. 
The  effort  to  partially  automate  the  free-space  bridge  dielectometer  is 
in  progress  at  the  present  time. 
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5 

FABRY-PEROT  RESONATOR  FOR  DIELECTRIC  MEASUREMENTS 


The  Fabry-Perot-type  resonator  is  very  useful  for  measuring 
dielectric  properties  of  materials  and  offers  a  high  degree  of  accuracy 
and  relative  ease  of  use.  The  measurements  can  be  made  with  the  same 
flat  panel  samples  used  in  the  free-space  bridge  dielectometer ,  making 
the  resonator  useful  for  providing  accurate  room-temperature  reference 
data  for  the  high-temperature  measurements  described  in  the  previous 
section. 

The  resonator  is  configured  as  a  semiconfocal  resonator:  one 

plane  mirror  and  one  spherical  mirror.  This  configuration  is  less 
sensitive  to  a  slight  misalignment  of  the  two  mirrors  than  for  the  case 
with  two  plane  mirrors.  Both  mirrors  are  machined  from  brass.  There 
is  a  single  iris  in  the  center  of  the  plane  mirror  to  couple  into  the 
resonator.  The  spherical  mirror  can  be  moved  along  the  symmetry  axis 
of  the  resonator  using  a  precision  micrometer  drive  that  reads  out  to 
0.001  millimeter.  Figure  5-1  is  a  schematic  diagram  of  the  resonator 
showing  dimensions  that  pertain  to  the  measurement  of  dielectric 
constant.  Referring  to  the  diagram  R0,  the  radius  of  curvature  of 
the  spherical  mirror  is  5.00  inches.  L  is  the  length  of  the  cavity  at 
resonance  for  a  particular  mode  without  the  sample  in  place.  Ah  is  the 
mirror  displacement  required  to  restore  resonance  at  the  same  mode 
after  the  sample  is  placed  in  the  resonator. 

With  the  sample  placed  against  the  flat  mirror,  as  shown  in 
Figure  5-1,  the  relative  dielectric  constant  tr  is  obtained  from  the 
index  of  refraction,  which  is  given  by  (Reference  5-1) 


tan  (nBd)  _  +  tan  8(d  +  AL)  . ,  . 

nBd  '  Rd  °  U 

where  n  =  /e^  and  B  =  2n/X0.  The  equation  is  solved  for  n  by 

numerical  iteration. 
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FIGURE  5-1.  Schematic  Diagram  of  the 
Semiconfocal  Resonator. 


The  loss  tangent  of  a  sample  is  determined  by  Pleasuring  the  qual¬ 
ity  factor  Q0  of  the  empty  resonator  and  Q  of  the  resonator  with  the 
sample  in  place.  Quite  generally,  the  Q  of  a  cavity  depends  on  the 
strength  of  the  coupling  into  the  cavity.  The  Q  of  a  cavity  is  defined 


stored 


energy  lost 


per  cycle 


As  the  coupling  is  increased,  the  oscillation  in  the  cavity  is  more 
rapidly  damped  and  therefore  Q  is  decreased.  Placing  the  test  sample 
directly  against  the  coupling  iris  in  the  flat  mirror,  as  shown  in 
Figure  5-1,  changes  the  coupling  of  the  resonator.  Since  we  are  inter¬ 
ested  in  measuring  the  change  in  Q  of  the  resonator,  which  is  produced 
solely  by  loss  in  the  sample,  the  sample  must  be  spaced  some  distance 
away  from  the  flat  mirror  to  avoid  affecting  the  coupling.  At  35  giga¬ 
hertz  we  have  obtained  adequate  results  with  spacing  as  little  as  0.065 
inch  off  of  the  flat  mirror. 
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The  Q  of  the  cavity  can  be  measured  by  noting  the  width  between 
half-power  points  on  the  resonance  line  as  the  micrometer  is  turned  and 
using 


Q 


L_ 

AL 


where  L  is  the  length  of  the  cavity  and  AL  is  the  3-decibel  displace¬ 
ment.  AL  can  be  measured  with  precision  up  to  four  significant  digits. 
Much  more  accurate  is  the  frequency  change  method  where  Q  is  given  by 


Q 


f 

o 


Af 


The  3-decibel  linewidth  Af  can  be  read  with  precision  up  to  nine 
significant  digits,  depending  on  the  stability  of  the  RF  source. 

The  loss  tangent  of  a  sample  is  measured  by  the  following  proce¬ 
dures  : 

•  Obtain  resonance  in  the  empty  resonator  by  moving  the  spherical 
mirror  and  note  the  length  L,  of  the  cavity  at  resonance. 

•  Measure  Q0  of  the  empty  resonator  by  frequency  change. 

•  Insert  sample  into  the  resonator  using  spacers  to  space  the 
sample  away  from  and  parallel  to  the  plane  mirror. 

•  Return  resonator  to  resonance  at  the  same  TEMQOq  mode 
obtained  in  the  empty  cavity  by  reducing  the  length  of  the 
cavity  using  the  micrometer. 

•  Measure  Q  of  the  cavity  with  the  sample  by  frequency  change. 

For  samples  of  thickness  greater  than  \/2,  the  loss  tangent  is  given 
approximately  by  the  following  expression  (Reference  6-2): 


Zi/Q  -  Lj/Q 

tan  5  =  — - -  (  5-2 ) 

(A  ,/A  )ZK  d  +  1 
1  o  r 


where 
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d  -  AL 


(5-3) 


(M2  ,  _L±. 

\  A  /  1  +  n 


cot2  (ngd) 

2  cot2  (nSd) 


(5-4) 


If  d  =  A/2,  the  error  in  Equation  5-2  is  approximately  10%,  which  is 
comparable  to  the  typical  experimental  error  at  these  wavelengths. 

Equation  5-2  is  derived  by  volume  integrations  of  the  energy 
density  in  the  resonator,  so  (to  the  extent  that  curvature  of  the  wave 
front  incident  on  the  sample  can  be  ignored)  the  spacing  of  the  sample 
from  the  plane  mirror  is  unimportant.  In  actuality,  plane  waves  exist 
only  at  the  surface  of  the  plane  mirror;  thus,  for  a  sample  of  appreci¬ 
able  thickness  or  large  separation  from  the  plane  mirror.  Equations  5-1 
and  5-2  must  be  corrected  for  wave-front  curvature.  A  number  of  dif¬ 
ferent  approaches  to  calculating  the  corrections  have  been  worked  out 
including  the  methods  of  Cullen  and  Yu  (Reference  5-3),  Cook  and  Jones 
(Reference  5-4),  and  Lynch  (Reference  5-5).  At  the  present  time,  these 
corrections  are  smaller  than  the  typical  measurement  errors  that  are 
encountered  using  the  NWC  apparatus;  thus,  for  the  data  herein 
reported,  these  corrections  are  not  employed. 

Figure  5-2  shows  a  block  diagram  of  the  NWC  laboratory  setup  for 
dielectric  measurements  using  the  Fabry-Perot  resonator.  The  RF  source 
consists  of  a  microwave  sweeper  that  is  phase  -  locked  to  the  internal 
time  base  of  a  phase-locking  frequency  counter  (EIP  Model  578).  This 
arrangement  can  be  used  much  like  a  microwave/millimeter  wave  synthe¬ 
sizer,  since  the  output  frequency  of  the  sweeper  can  be  controlled  by 
front  panel  numerical  input  to  the  frequency  counter  with  a  resolution 
of  10  kilohertz.  Directional  couplers  are  used  to  sample  equally  the 
power  incident  on  and  power  reflected  from  the  resonator  coupling  iris 
at  the  end  of  the  waveguide  feed.  Matched  crystal  detectors  are  used 
to  measure  the  ratio  of  reflected/incident  power  that  is  displayed  on  a 
scalar  network  analyzer.  By  ratioing  the  reflected  and  incident  power 
levels,  measurements  are  made  less  sensitive  to  amplitude  fluctuations 
or  amplitude  modulation  noise  from  the  sweeper. 

The  resonator  itself  is  constructed  to  accept  interchangeable 
feeds  to  cover  multiple  frequency  bands,  currently  at  26.5  to  40  giga¬ 
hertz  and  60  to  90  gigahertz.  At  35  gigahertz,  the  resonator  exhibits 
a  maximum  loaded  Q  of  20  000.  At  the  present  time,  no  measurements 
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FIGURE  5-2.  Block  Diagram  of  Setup  for  Dielectric  Measurements. 

have  been  made  in  t lie  higher  frequency  band.  Figure  5-3  is  a  photo¬ 
graph  of  the  resonator  with  a  dielectric  sample  in  place  and  spacers 
used  for  loss  measurements.  Figure  5-4  shows  the  interchangeable  wave¬ 
guide  feeds  that  are  used  with  the  resonator. 
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Measurements  of  dielectric  constant  and  loss  tangent  were  made  at 
35  gigahertz  in  the  Fabry-Perot  resonator  for  comparison  with  data 
obtained  in  the  free-space  bridge  dielectometer  (Table  5-1).  The 
results  from  the  resonator  measurements  are  presented  in  the  table. 


TABLE  5-1.  Resonator  Measurements. 


Material 

e 

r 

tan  5 

Duroid  5880 

— 

— 

Duroid  5870 

2.22 

0.0012 

1-Cerez 

1.98 

0.0056 

2-Cerez 

1.94 

0.0021 

Pyroceram  9606 

5.63 

0.00026 

At  present,  the  sensitivity  of  the  measurements  made  in  the 
Fabry-Perot  resonator  apparatus  is  limited  by  the  stability  of  the  RF 
source.  By  employing  a  phase-lock  loop  with  wider  loop  bandwidth  and  a 
more  stable  reference,  the  RF  frequency  could  be  controlled  to  a  reso¬ 
lution  of  100  hertz  or  better.  Addition  of  an  automatic  gain  control 
for  the  RF  source  could  also  enhance  the  measurement  precision. 
Synchronous-detection  and/or  signal-averaging  techniques  would  allow 
the  possibility  of  looser  coupling  to  the  cavity,  thereby  increasing 
the  Q  factor.  Much  work  remains  to  be  done  in  terms  of  acquiring  data 
on  materials  in  the  35-gigahertz  range  and  above.  This  effort  shall 
proceed  as  more  materials  become  available  for  testing. 
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6 

COMPUTER-CONTROLLED  MOVING  VANE  DIELECTRIC 
CONSTANT  MEASUREMENT  SYSTEM 


The  complex  propagation  constant  y  and  characteristic  wave  imped¬ 
ance  Z  of  a  material  are  determined  by  the  complex  relative  dielectric 
permittivity  e*  =  e( 1-j  tan  6e)  and  magnetic  permeability  p*  = 
u(l-j  tan  Sm) .  F0r  free-space  propagation,  the  relations  are 


Y 


(6-la) 


Z  =  Z 

o 


(6-lb) 


For  a  material  of  thickness  i,  backed  by  a  conducting  surface,  the 
complex  reflection  coefficient  r(  £)  is  given  by 


rU) 


r(  “)  -  e 


-2  y£, 


i  -  r(-) 


-2y  l 


( 6-2a) 


where 
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The  quantity  T(")  represents  the  limit  yfc  >  1  and  corresponds  to  the 
empty  space/material  interface  reflection  coefficient.  In  guided-wave 
structures,  the  parameters  y  and  Z  are  determined  by  e*  and  p*  as 
well  as  the  frequency  and  waveguide  dimensions  (see  Appendix  B). 

The  moving  vane  dielectometer  is  illustrated  in  Figures  6-l(a)  and 
(b).  The  moving  vane  is  a  solid  metal  bar  that  completely  fills  the 
waveguide.  A  slot  is  cut  in  this  bar  to  accommodate  the  sample  to  be 
measured.  The  width  of  this  slot  is  made  sufficiently  small  such  that 
it  represents  a  waveguide  well  below  cutoff  when  loaded  with  the 
dielectric  sample.  The  sample  is  centered  in  the  waveguide  and  remains 
stationary  during  a  measurement.  The  vane  is  driven  by  a  stepper-motor 
assembly.  Figure  6— 1(b)  illustrates  three  positions  of  the  vane  rela¬ 
tive  to  the  sample.  For  the  upper  drawing  (sample  recessed),  the  vane 
approximates  a  moving  short.  In  the  lower  drawing,  the  vane  is 
recessed  and  the  sample  inhomogeneously  fills  the  waveguide.  The 
"effective"  sample  length  is  determined  by  the  vane  position  and  is 
varied  continuously  during  a  measurement.  A  detailed  analysis  of  the 
inhomogeneously  filled  waveguide  is  presented  in  Appendix  C. 
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FIGURE  6-1. 
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Sample  and  Moving  Vane  Assembly. 
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FIGURE  6-1-  Contd. 


A  block  diagram  of  the  measurement  system  is  shown  in  Figure  6-2. 
Measurement  errors  arise  from  the  non-square-law  power  dependence  of 
the  mixer  crystals  (see  Appendix  D)  and  their  power-level-dependent 
direct  current  offsets.  The  former  is  eliminated  by  inserting  20-  to 
30-decibel  attenuation  (step  attenuator)  in  the  signal  arm.  To  elimi¬ 
nate  direct  current  offsets,  four  independent  measurements  are  made 
versus  sample  position  at  the  relative  phases  (precision  phase  shifts) 
0,  90,  180,  and  270  degrees.  Subtracting  data  pairs  (0  to  180  degrees 
and  90  to  270  degrees)  provides  the  corrected  in-phase  and  quadrature 
components  of  the  reflection  coefficient.  This  procedure  is  performed 
for  each  vane  position  in  a  measurement  set. 

The  moving  vane  system  has  been  automated  and  is  controlled  by  an 
HP- 1000/F  computer.  A  special  digital  interface  has  been  designed  that 
controls  stepper-motor  acceleration,  synchronizes  stepper-motor  posi¬ 
tion  and  data  acquisition,  and  repositions  the  vane  to  precise  initial 
positions.  For  each  of  the  reference  phase  values,  the  vane  is  cycled 
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FIGURE  6-2.  Block  Diagram  of  the  Moving  Vane 
Dielectric  Measurement  System. 


from  minimum  to  maximum  to  minimum  vane  insertion  positions.  Data  for 
each  cycle  are  stored  as  a  disk  file  and  the  subtractions  are  performed 
with  software.  The  total  sample  measurement  time  is  under  4  minutes. 

A  series  of  measurements  on  different  representative  materials  has 
been  performed.  The  measurement  conditions  were  F  =  9.7  gigahertz  and 
sample  (slot)  widths  of  0.2  inch  in  rectangular  X-band  waveguide. 
Results  for  a  6-guide-wavelength-long  Teflon  sample  gives  e  =  1.966. 
Individual  half-wavelength  segments  were  analyzed  and  averaged  over  the 
sample  length  giving  a  root-mean-square  deviation  of  0.64%. 

Data  for  a  Bakelite  sample  were  obtained  and  are  plotted  as 
sin(phase)  and  phase  versus  position  in  Figures  6-3(a)  and  (b),  respec¬ 
tively.  The  "empty"  and  "inhomogeneously  filled"  regions  of  relative 
vane  position  are  readily  apparent  from  the  period  differences.  The 
distortion  in  the  "filled"  region  is  caused  by  the  standing  wave  estab¬ 
lished  between  the  empty/filled  waveguide  interface  and  the  vane  posi¬ 
tion.  The  standing  wave  contribution  is  illustrated  in  Figure  6—  3(b) 
for  phase  versus  position.  In  the  absence  of  the  interface/vane  stand¬ 
ing  wave  contributions,  the  empty/filled  regions  would  have  constant 
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slopes  inversely  proportional  to  the  guide  wavelength.  The  result  for 
this  sample  is  e  =  3.054.  The  root-mean-square  standard  deviation  of 
2.2%  compared  to  the  result  of  0.6%  for  Teflon  indicates  the 
inhomogeneity  of  the  dielectric  constant. 

Figure  6-4  is  the  sin(phase)-versus-position  variation  for  a  metal 
foam  artificial  dielectric  microwave  lens  material.  A  close  inspection 
of  this  figure  indicates  substantial  volume  inhomogeneity  in  the 
dielectric  constant.  The  average  value,  e  =  2.96,  exhibits  a  root- 
mean-square  deviation  of  3.8%  from  the  mean. 


FIGURE  6-4.  Sin(Phase)  Variation  for  Metal  Foam  Artificial 
Dielectric. 


Figures  6-5(a)  and  (b)  are  initial  results  for  a  special  composite 
magnetic  absorber.  The  results  in  Figure  6-5(a)  show  the  magnitude  of 
the  reflection  coefficient  versus  sample  length  for  a  rectangular  sam¬ 
ple.  The  results  in  Figure  6-5(b)  show  the  same  material  tapered  from 
zero  to  full  height  over  a  length  of  approximately  3  guide  wave¬ 
lengths.  Effective  elimination  of  the  interface  reflection  will  facil¬ 
itate  determination  of  absorption  coefficients  in  this  class  of  highly 
absorbing  materials. 
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Appendix  A 

COMPUTER  PROGRAM  FOR  NUMERICAL  INVERSION 
OF  FLAT  PANEL  TRANSMISSION  COEFFICIENT 
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LISTING  OF:  &DIELC 

MADE  BY  GAGNON . GENERAL  AT  9;02  AM  TUE  .  ,  3  DEC. 


1985 


0000  1 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00015 
00016 
00017 
00018 
00019 
00020 
00021 
00022 
00023 
00024 
00025 
00026 
00027 
00028 
00029 
00030 
00031 
00032 
00033 
00034 
00035 
00036 
00037 
00038 
00039 
00040 
00041 
00042 
00043 
00044 
00045 
00046 
00047 
00048 
00049 
00050 
00051 
00052 
00053 


FTN7X 


01 


05 


DIELC  version  of  11/25/85 

This  program  is  used  to  determine  the  complex  permittivity 
of  a  single-layer  flat  panel  of  dielectric  material  from 
a  measurement  of  the  phase  and  amplitude  of  transmitted 
power  through  the  sample  for  any  angle  of  Incidence.  The 
program  employs  Mu 1 1 er ' sme t hod  for  determining  the  zeroes 
of  a  function  of  a  complex  variable. 
********************************************************* 


program  dielc 

common  trans,  dbeta,  theta,  mpol 
complex  xg,  xz,  trans 
real  eps,lamda 

write!!,'!" ENTER  ESTIMATED  DIELECTRIC  CONSTANT")' 

read ( 1 , * )esubr 

xg- cmplx ( esubr ,  0.) 

nd-  7 

eps-  1 . e-8 
nit-  100 

wr ite! 1 ,'( "ENTER  FREQUENCY  IN  GHZ" ) ' ) 
read ( 1  ,  *  )  f  req 


f  req- f  req* 1 . e9 
vrite( 1 , * ( "ENTER  1 
write! 1 , ' ( "ENTER  2 
read! 1 , *  Jmpol 
i f ( mpo 1 . ne . 1  . and . 


FOR  PARALLEL  POLARIZATION")1) 

FOR  PERPENDICULAR  POLARIZATION")*) 


mpo 1 . ne . 2 )  go  to  1 
write! 1 ,  *( "ENTER  SAMPLE  THICKNESS  IN  INCHES")') 
read ( 1  ,  *  )  d 

wr lte! 1 ,*( "ENTER  INCIDENCE  ANGLE  IN  DEGREES" )* ) 
read! 1 , * ) theta 

wr ite! 1 ,'( "ENTER  MEASURED  PHASE  SHIFT  IN  INCHES"  )  1  ) 
read! 1 ,* )deltaphi 

wri te( 1 ,  '< "ENTER  TRANSMISSION  MAGNITUDE")*) 

read (  1  ,  *  ) tmag 

pi -atan ( 1 . ) *4 . 

lamda-2 .998el0/(2.54*freq) 

dbeta-d*2 . *pi /lamda 

theta- t he ta*pi/180. 

phi-deltaphi+d*cos(theta) 

phl-2.*pi*phi/l amda 

trans  is  the  measured  transmission  coefficient. 

trans-tmag*cexp( (0. , -1 . )*phi ) 

call  muller(xg,  nd, eps, nit,  xz,iwhy) 

write!!, ' C  " ) * ) 

wr ite! 1 ,* ) *  EPSILON- *  ,xz 

tandelta--aimag(xz)/real(xz) 

write! 1 ,*) *  LOSS  TANGENT- *  ,tandelta 

l f ! 1 why . ne . 3 )  go  to  5 

wr i te ( 1  ,  *  ( " NO  CONVERGENCE ! " ) ’  ) 

conti nue 

end 

subroutine  mul ler(xg,nd,eps,nit,xz, iwhy ) 
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0005*  c 

0005  5  c  We  follow  Atkinson,  ''Intro,  to  numer leal  analysis''  ,  sec.  2.10 

00056  c  Recall  that  one  is  allowed  to  permute  indices  at  will,  sort  of, 

00057  c  in  Wewton's  divided-differences . 

00058  c  fen  needed 

00059  c  no  deflation  (only  i  zero  gotten) 

00060  c  input  : 

00061  c 

00062  c  xg, 

00063  c  nd  : 

00064  c 

00065  c  eps  s 

00066  c  nit, 

00067  c 

00068  c 

00069  c  output: 

00070  c  xz ,  the  answer  (if  iwhy  not  3) 

00071  c  iwhy,  the  reason  it  guit.  See  above. 

00072  c 

00073  common  trans,  dbeta,  theta,  mpol 

00074  complex  xg,xz,  x(0:3),h,  f21,f20,f0i,  f(0:3),  f012, 

00075  *  w,  dent, den 

00076  c 

00077  c  x(3)  is  the  new  estimate,  based  on  x0,xl,x2,  with 

00078  c  x(2)  the  old  estimate 

00079  real  eps,  epsx 

00080  integer  nd ,  nit,  ic,  iwhy 

00081  complex  fen 

00082  external  fen 

00083  logical  convl 

00084  c 

00085  epsx  •  l./lo.**nd 

00086  ic  -  0 

00087  h  •  (0.01,  0. ) 

00088  if  (cabs(xg)  .gt.  0.1)  h-(0.1,  0.)*xg 

00089  x(0).  xg-h 

00090  x(l)-  xg+h 

00091  x  (  2  )  -  xg 

00092  f(0)  •  f  cn ( x ( 0 ) ) 

00093  f(l)  -  f cn ( x ( 1 ) ) 

00094  f ( 2 )  -  fcn(x(2)) 

00095  f 21  -  (f ( 1  )-f (0) )/(x( 1 )-x( 2  )  ) 

00096  10  continue 

00097  fOl  -  f  2 1 

00098  c  - Recall  that  f  [  xO , xl ] - f ( xl , xO ]  ,  in  Newton's  world. 

00099  f  2 1  -  <f ( 1 )-f ( 2) )/(x( 1 )-x( 2 )  ) 

00100  f 20  -  (f ( 2 ) -f (0) )/(x(2)-x(0)  ) 

00101  w  -  f  2 1  +f  20  -fOl 

00102  f  0 1 2  -  ( f  2 1  -  fOl )/(x(2)-x(0) ) 

00103  dent  •  csgrt(w»w  -  4 . 0* f ( 2 ) »f 0 1 2 ) 

00104  if  ( cabs ( w+den t )  ge.  cabs ( w-dent ) )  then 

00105  den-  w+dent 

00106  else 

00107  den-  w-dent 

00108  endif 

00109  x ( 3 )  -  x ( 2 )  -  (2 ,0*f ( 2) )/den 


the  guess  at  x  with  fcn(x)-  0. 

convergence,  If  nd  digits  agreement  in  two  successive 
estimates  of  x  (iwhy  set  to  1) 

convergence,  if  x  found  within  eps  (iwhy  set  to  2), 
the  maximum  number  of  iterations  allowed  (Iwhy  set  to  3, 
if  nit  iterations  occured  without  convergence) 
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00110 
00111 
00112 
0011  3 
00114 
00115 
00116 
00117 
00118 
001  19 
00120 
00121  20 
00122 
00123 
00124 
00125 
00126 
00127 
00128 
00129 
00130 
00131 
00132  c 
00133 
00134 
00135 
00136 
00137  c 
00138  c 
00139 
00140 
00141 
00142 
00143 
00144 
00145  c 
00146 
00147 
00148 
00149  c 
00150  c 
00151 
00152 
00153 


f (3)  •  f  cn ( x ( 3 )  ) 

1C  -  1C  +  1 

convl -  cabs(x(3)-x(2) ) . le. epsx*amaxl ( 0.01,cabs(x(3))) 
if  ( convl . or . cabs ( f ( 3 ) ) . le . eps  .or.  ic.ge.nit)  go  to  20 
x  (  0  )  -  x  (  1  ) 
x  (  1 )  -  x( 2  ) 
x  (  2  )  -  x  (  3  ) 
flO)  -  f(l) 
f(l)  -  f ( 2  ) 
f(2)  -  f  (  3  ) 
go  to  10 
continue 
xz  •  x( 3 ) 
if  ( convl )  then 
iwhy- 1 

elseif  ( cabs ( f ( 3 ) ) . le . eps  )  then 
iwhy- 2 
else 
iwhy- 3 
endi  f 
return 
end 

******************************************************* 

complex  function  fcn(x) 

common  trans,  dbeta,  theta,  mpol 

complex  x,t, trans 

complex  radsq,rad,rpe,rpa,r,psi ,psiexp 

calculate  a  value  for  interface  reflection,  r,  . . 

r ad sq*x-sin( theta )**2 
rad-csqr t ( radsq ) 

rpa-(x*cos(theta)-rad)/(x*cos(theta)+rad) 
rpe-(cos(theta)-rad)/(cos(theta)+rad) 
if (mpol . eq . 1 )  r-rpa 
if (mpol . eq . 2 )  r-rpe 

- then  calculate  a  value  for  transmission,  t. 

ps i -dbeta*rad 

psiexp-cexp( (0. , -1 . ) *psi ) 

t- (1 . -r*r )*psiexp/( 1 . -r *r *ps iexp**2 ) 

fen  to  be  zeroed  is  the  difference 

f cn- t rans-t 

return 

end 
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Appendix  B 

ERROR  MODEL  FOR  MULTIPLE  REFLECTIONS  IN  THE 
FREE-SPACE  BRIDGE  DIELECTOMETER 


In  Figure  B-l,  T  is  the  complex  transmission  coefficient  of  the 
system  with  the  sample  in  place  and  is  the  actual  measured  quantity. 
The  transmission  magnitude  t  and  transmission  phase  $  of  the  sample  are 
the  desired  quantities  for  determining  the  sample  dielectric  constant 
and  loss  tangent,  respectively.  R  is  the  reflection  coefficient  of  the 
lens  aperture  of  each  horn,  assuming  they  are  identical.  The  quantity 
r  represents  the  reflection  from  the  horn  feed  plus  any  small 
reflections  appearing  down  the  transmission  line,  and  y  is  an  unknown 
phase  constant.  The  three  unknowns,  R,  r,  and  y,  are  determined  by 
least-squares  fit  to  data  for  transmission  phase  or  magnitude  versus 
separation  of  the  two  horns  as  measured  from  the  micrometer-screw 
positioner.  For  the  purpose  of  error  correction,  the  two  different 
reflections  can  actually  be  lumped  together,  greatly  simplifying  the 
model. 


-  =  _B_  =  U  -  r~)  11  -  R  )  1  e  '  e  lv _ 

A  '(1  +rRe-l2V  _  (Rf,.-*2*,2,*,-'2* 

T  =  iTl  e-*'*’ 

sin  (2y  -  0)  -  t2  (R2  +  r2|  s,n  I2y  +  0)  -  2r  R  (1  +  t2)  sin  0 
tan  <X>  -  _ _ _ 

cos  (2y  -  0)  -  t2  (R2  +  r2|  cos  (2->  +  O)  +  2r  R  (1  -  t2)  cos  0 
ITl  *  I  (1  +  4  r  R  cos  2>  -  2t2  (R2  +  r2)  cos  0  -  4  r  R  t2  cos  (2i  -  20)]  ~1/2 


FIGURE  B-I.  Error  Model  for  Multiple  Reflections  in  the  Free-Space 
Bridge  Dielectometer . 
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Appendix  C 

MULTIPLE-DIELECTRIC-SECTION  INHOMOGENEOUSLY  FILLED  WAVEGUIDE 


Composite  dielectric  sections  in  rectangular  waveguides  have  been 
treated  by  Marcuvitz  (Reference  C-l)  and  Collin  (Reference  C-2).  For 
TEmo  modes,  the  electric  field  is  parallel  to  the  dielectric  sur¬ 
faces  and  the  propagation  constant  can  be  determined  from  a  transverse 
equivalent  transmission  line  circuit.  Solutions  for  two-  and  three- 
dielectric  sections  have  been  obtained  by  this  method  and  are  shown  in 
Figure  C-l.  Applied  to  the  analysis  of  actual  dielectric  sections  in 
waveguide,  the  perfect  filling  shown  in  Figure  C-l  is  an  approxima¬ 
tion.  Gaps  between  the  dielectric  surfaces  and  the  metallic  walls  and 
between  the  dielectric  sections  exist  in  laboratory  realizations  of 
these  structures.  Describing  these  "gaps"  as  additional  dielectric 
sections,  the  ideal  two-  and  three-section  problems  in  Figure  C-l 
become  five-  and  seven-section  problems,  respectively.  A  general  solu¬ 
tion  for  TEmo  modes  with  N  dielectric  sections  is  presented  here. 
The  interim  sections  are  represented  as  2  x  2  matrices  combined  through 
N-2  matrix  multiplications.  Whether  two,  three,  or  N  sections  are  con¬ 
sidered,  the  resulting  transcendental  equation  must  be  solved  numeri¬ 
cally.  The  general  solution  has  been  implemented  as  a  function  sub¬ 
routine  requiring  as  input  the  width  and  dielectric  constant  of  each 
section  and  the  number  of  sections.  Its  advantage  is  that  different 
dielectric  section  multiplicities  can  be  analyzed  with  a  single  proce¬ 
dure. 


FIGURE  C-l.  Configuration  of  the  Ideal  Two-  and  Three-Section  Inhomo- 
geneously  Filled  Waveguides. 
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The  N-section  dielectric-filled  waveguide  is  shown  in  Figure  C-2. 


FIGURE  C-2.  Configuration  of  the  N-Section  Dielectric- 
Filled  Waveguide. 


Restricted  to  TEmo  modes,  the  electric  field  in  the  nth  section  can 
be  written  as 


E  (n)  =  D  cos  k  x  +  D'  sin  k  x 
y  n  n  n  n 


( C- 1  a  ) 


where 


(C-lb) 


(C-Lc) 


The  quantity  en  is  the  relative  dielectric  constant  of  the  nth  sec¬ 
tion.  For  materials  with  nonzero  dielectric  loss,  en  =  e'n 
je''n.  Continuity  of  Ey  and  Hz  at  the  nth  interface  can  be 
represented  by 


M.(k  ,  x  )  •  6  =  M  (k  ,  x  )  •  i5 

I  n  ’  n  n  r  n+1  n  n+1 


(C-2) 
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+  rj> 

where  Dn  =  (Dn,  D'j,)1.  The  field  amplitude  in  the  1st  and 
nth  sections  are  related  by 


Mg(kl»  xl)  *  ^1  =  f (k 2 ,  xj)  •  M  *(k2,  x2)  1 


'VVr  Vz>  '  "I'Vr  W 


W  Vl>  •  ®» 


(C-3) 


Defining  the  matrices  M’(n)  and  S  as 


M'(n)  =  W  Vl>  *  M^(kn*  Xn> 


(C-4a) 


S  =  M 1 ( 2)  •  M * ( 3 )  ...  M'(N-l) 


(C-4b) 


the  field  amplitude  relations  become 


Vkl>  xl>  *  =  S  *  W  xN_j)  *  SN 


(C-4c) 


The  components  of  M'(n)  are 


M’(n)  = 


cos  k  (x  -  x  -  k  1  sin  k  (x  -  x  ,) 

n  n  n-1  n  n  n  n-1 


<•  sin  k  (x  -  x  cos  k  (x 

n  n  n  n-1  n 


x  ,  ) 
n  n  n-1 


(C-5) 


The  electric  fields  are  null  at  the  metallic  boundaries  xq  and  x 
and  give  the  conditions  ^ 
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Mr(kN’  ’Sj-l*  *  °N 


'Sin  kN(XN  “  XN-1} 


kN  cos  kN(xN  -Vl) 


D(N) 


(C-6a) 


and 


M £(k  1 »  M)  •  Dn 


-  sin  kj(xj  -  x 0 ) 

-  k  j  cos  k  j (x  j  -  Xg) 


D(l) 


(C-6b) 


where  D(N)  and  D(l)  are  scalars.  Substituting  into  Equation  C-4c  gives 
two  equations  in  D(N)  and  D(l).  A  nontrivial  solution  requires 


S2  1 

FTiT  sln  ki(xi 


X0)  sin  kjjtXjj  -  xN_j) 


+  S12  cos  k ! (x  j  -  x0)  cos  kN(xN  -  x^) 


-  sin  kjCxj  -  x 0 )  cos  kN(xN  -  x^) 


11 


+  ~  cos  kl(xl  "  x0>  sin  kN(xN  ~ 
N 


=  0 


(  C— 7 ) 


Equation  C-7  represents  the  general  solution  for  the  N-section  problem. 

The  general  solution  can  be  represented  more  compactly.  Defining 
Tn  =  tan  kn(xn  -  xn_i),  the  general  transcendental  equation  be¬ 
comes 


NWC  TP  6643 


titn  ti  tn 

S21  kikN  S22  kj  +  Sl1  kN  +  S12 


(C-8) 


Since 


M'(n)  = 


(C-9) 


for  numerical  calculations,  the  matrix  of  a  section  (Equation  C-9)  is 
determined  and  multiplied  by  S,  producing  the  new  value  of  S 
(Equation  C-4).  This  procedure  is  repeated  for  each  of  the  N-2 
interior  sections.  The  results  are  substituted  into  the  "function” 
(Equation  C-8).  Zeros  of  Equation  C-8  and  thus  solutions  are 
determined  using  standard  numerical  techniques.  In  practice,  the 
functional  form  of  Equation  C-7  is  used,  avoiding  the  numerical 
problems  associated  with  poles  in  the  tangent  functions. 

Applying  the  general  result  of  Equations  C-8, C-4,  and  C-9  to  the 
two-,  three-,  and  four-section  problems  gives 


Two  Sections 


Three  Sections 


(k2  T; 


7— "  +  ; — —  +  =  0 


Four  Sections 


+  k3T3)  +  ( l-3  t2t3^  +  ^  t2t,\ 
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The  symmetrical  three-section  solution  ej  =  e3,  x3  =  a,  and  x2  -  Xj  =  d 
gives 


The  antisymmetric  solution  (tan,  tan)  differs  by  a  sign  from  the  result 
given  by  Collin  (Reference  C-3). 

As  noted  at  the  beginning  of  this  section,  the  general  solution 
(Equation  C-8)  is  readily  implemented  numerically  as  a  single 
procedure,  which  makes  explicit  expressions  for  the  cases  of  different 
multiplicities  unnecessary.  It  should  also  be  noted  that  as  the 
frequency  is  increased  above  cutoff,  at  some  point  the  angle  of 
incidence  for  a  particular  layer  will  exceed  the  critical  angle  in 
propagating  from  a  higher  dielectric  material  to  a  lower  one 
(Reference  C-4).  This  leads  to  total  internal  reflection.  In  such  a 
case,  an  evanescent  wave  exists  in  the  lower  dielectric  constant 
material,  and  the  tangent  function  in  Equation  C-8  becomes  a  hyperbolic 
tangent.  Further  complexity  results  in  lossy  materials  caused  by  the 
fact  that  real  and  imaginary  characteristic  equations  must  be  solved 
simultaneously. 
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Appendix  D 

BALANCED  MIXER  MAGNITUDE  ERRORS:  DEVIATIONS  FROM  SQUARE-LAW  RESPONSE 

Analysis  of  data  obtained  with  a  moving  vane  dielectometer  system 
exhibits  magnitude  errors  that  cannot  be  attributed  to  the  finite 
reflection  coefficients  of  the  microwave  system  components  used.  The 
spectral  purity  of  the  power  source  was  measured  and  is  not  the  source 
of  these  errors.  Detector  crystal  square-law  deviations  become  the 
major  error  source  when  used  with  a  nearly  ideal  circulator. 

The  moving  vane  dielectometer  is  computer  controlled.  Data  are 
obtained  as  follows:  For  the  phase  shifter  in  the  mixer  reference  arm, 
set  first  at  0  degree  then  at  90  degrees;  the  in-phase  and  quadrature 
variations  versus  vane  position  are  determined  and  stored  as  disk 
files.  The  position  dependence  of  the  magnitude  and  phase  of  the 
reflection  coefficient  are  computed.  Data  obtained  for  a  moving  short 
at  9.7  gigahertz  are  shown  in  Figure  D-l.  In  the  absence  of  system 
errors,  the  magnitude  should  be  a  constant. 


FIGURE  D-l.  Measured  Magnitude  Variation  of  the  Voltage  Reflection 
Coefficient  for  a  Moving  Short  in  the  Waveguide;  Nonideal  Bahavior. 
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The  microwave  circuit  (Figure  6-2)  has  been  analyzed.  Nonideal 
circulator  response  is  represented  by  (T22,  T13  *  0)  in  the  response 
matrix  T.  The  output  of  the  mixer  is  represented  by 


where  Eg  and  Er  represent  the  electric  fields  (complex)  from  the 
"signal";  Pi  and  n2,  the  niixer  crystal  square  law  deviation;  and  e, 
crystal  unbalance. 

The  "apparent"  magnitude  of  the  moving  short  reflection  coeffi¬ 
cient  has  been  calculated  numerically  and  is  plotted  versus  the  phase 
(relative  position)  of  the  reflection  coefficient  with  no  system 
errors . 

Figure  D-2  shows  the  behavior  of  parameters  that  represent  the 
general  form  of  the  experimentally  measured  variation  (Figure  D-l)  and 
corresponds  to 


T22  =  0.02,-45°,  T13  =  0 
(hi,  n2)  =  -0.2 


e  =  0.01 


From  data  obtained  with  a  moving  short,  the  parameters  character¬ 
izing  the  square-law  deviations  of  the  mixer  crystals  can  be  determined 
using  conventional  numerical  curve-fitting  techniques.  Alternately, 
the  measured  moving  short  response,  tabulated  as  magnitude  versus 
phase,  can  be  used  directly  to  correct  system  errors.  The  power  level 
dependence  of  square-law  deviations  requires  the  use  of  identical 
average  powers  at  the  mixer.  Clearly  preferable  to  this  approach  is 
a  power  decrease,  so  that  square-law  deviations  are  both  minimized  and 
become  insensitive  to  small  power  changes. 
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